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SECTION  I 
INTRODUCTION 


A  detailed  analysis  of  the  requirements  for  an  electronic  shadowgraph  system  in  the 
Aeroballistics  Research  Facility,  located  at  the  Eglin  Air  Force  Base,  Ft.  Walton  Beach,  Florida, 
was  performed  in  1985  and  presented  at  the  International  Congress  on  Instrumentation  in 
Aerospace  Simulation  Facilities  (ICIASF)  meeting  (Reference  1)  at  Stanford  University  in 
August,  1985.  Briefly,  the  facility  is  680  feet  long  and  12  feet  wide  for  the  first  one  third  of  the 
range,  and  16  feet  wide  for  the  last  two  thirds.  Projectiles  and  models  of  various  shapes  are 
fired  down  the  range  and  photographed  at  50  stations  along  the  range  to  determine  their  stability, 
aerodynamic  properties  and  other  flight  characteristics.  Each  imaging  station  along  the  range, 
shown  in  Figure  1,  captures  two  orthogonal  shadowgraph  images  of  the  projectile  as  it  passes. 
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The  objective  of  the  analysis  was  to  provide  system  requirements  and  a  conceptual  design 
of  a  state-of-the-art  high-speed  electronic  imaging  system  to  improve  the  data  acquisition  and 
analysis  of  projectiles  in  flight.  The  projectiles  vary  in  shapes  and  may  travel  at  speeds  up  to 
MACH  6.  The  required  resolution  was  determined  by  measuring  the  transfer  function  and 
optical  illumination  parameters  of  the  existing  system. 

From  the  transfer  function,  the  limiting  resolution  was  determined  to  be  0.72  cycles  per 
millimeter  at  the  reflective  screen.  Illumination  parameters  were  determined  to  ascertain  the 
usefulness  of  the  available  spark-gap  light  source  for  large  area  arrays.  Intensity  and  spectral 
characteristics  of  the  source  were  measured. 

Based  on  the  analysis,  a  prototype  was  developed  to  verify  the  analysis  and  to  determine  if 
projectile  points  of  interest  can  be  located  with  the  desired  accuracy.  The  prototype  system 
described  in  Sections  II  and  HI  acquires  images  at  only  one  station.  Static  projectile  images 
were  successfully  acquired  using  modified  Westinghouse  vidicon  cameras  and  a  VICOM  image 
processing  system.  However,  the  vidicon  cameras  were  not  sensitive  enough  to  image  with  the 
Eglin  spark-gap. 

A  re-evaluation  of  the  project,  based  on  the  failure  of  the  vidicon  cameras  and  associated 
high  cost  of  comparable  charge-coupled  device  (CCD)  cameras,  resulted  in  evaluating  various 
algorithms  for  making  projectile  measurements  which  would  reduce  the  resolution  requirements. 
Two  of  those  algorithms  are  described  in  Section  V. 

Section  VI  describes  the  method,  and  derives  the  equations,  for  locating  a  projectile  in  3- 
space  from  the  centroid  and  second  moment  of  inertia  data  taken  from  two  orthogonal 
shadowgraphs  with  an  electronic  camera. 

As  a  result  of  the  promise  of  the  resolution  reducing  algorithms,  a  proof-of-principle 
imaging  system  was  assembled  and  installed  at  the  Eglin  Aeroballistics  Research  Facility. 

Images  on  various  shots  were  successfully  acquired  using  the  Eglin  spark-gap,  and  a  few 
measurements  were  made  to  compare  accuracy  with  the  film  cameras.  Section  VII  describes  that 
system  and  the  results  of  comparisons  with  film  measurements. 


SECTION  II 


EARLY  PROTOTYPE 


The  early  prototype  system  consisted  of  two  Westinghouse  vidicon  cameras,  an  image 
processing  system,  a  spark-gap  light  source,  an  infrared  linear  array  motion  detector,  and  a  retro- 
reflective  screen  on  which  a  shadowgraph  is  formed. 

1.  CAMERAS 

A  large  area  array  detector  of  2000  elements  high  and  1000  elements  wide  will  almost 
satisfy  the  resolution  requirements  for  the  Eglin  range  (Reference  1).  The  cameras  initially 
chosen  were  Westinghouse  model  ETV-2000  slow  scan  vidicons,  modified  to  output  a  2048v  x 
1024h  image.  When  they  were  selected,  they  were  the  only  cameras  known  that  would 
approximate  the  resolution  requirements.  The  cameras  are  operated  in  an  "integrate  mode"  (open 
shutter)  in  total  darkness.  They  are  erased  just  prior  to  firing  the  projectile,  and  the  image  is 
acquired  when  the  infrared  motion  detector  senses  the  projectile  and  triggers  the  spark-gap  light 
source.  As  soon  as  the  light  source  is  triggered,  the  image  on  the  vidicon  is  scanned  in  about  0.5 
seconds  and  sent  as  single  frame  video  to  be  digitized  by  the  image  processor. 

2.  IMAGE  PROCESSOR 

The  image  processor  is  a  VICOM  Systems  model  1800,  which  is  a  microprocessor  based 
system.  The  system  has  image  memory  space  for  two  2048v  x  1024h  x  16  bit  images,  for  a  total 
of  4  mb.  The  image  processor  has  a  50  mb  hard  disk,  printer,  console,  and  a  color  monitor 
capable  of  displaying  a  512v  x  5 12h  pixel  color  image.  Image  processing  software  is  also 
included,  which  allows  a  variety  of  processing  operations  such  as  filtering,  image  arithmetic, 
morphological,  and  graphical,  to  be  performed  on  images  either  interactively  or  as  a  subroutine 
called  from  other  software. 

3.  SPARK-GAP  LIGHT  SOURCE 

An  open  air  spark-gap,  Hi-Voltage  model  SS55P,  with  a  capacitor  discharge  circuit 
charged  to  4800  vdc  is  employed  as  a  strobed  light  source.  The  spark  is  discharged  across  a 
2.2  mm  air  path,  which  is  small  enough  to  be  considered  as  a  point  light  source.  The  spark-gap 
is  positioned  adjacent  to  each  camera  and  casts  a  well-defined  shadow  on  the  reflective  screen. 
The  spectral  characteristics  were  reported  in  the  1985  ICIASF  paper  (Reference  1).  The  output 
of  the  spark-gap  has  a  duration  of  about  300  nanoseconds.  This  duration  is  short  enough  to 
"stop"  a  projectile  travelling  at  Mach  6. 

4.  INFRARED  ARRAY  TRIGGER 

A  narrow  infrared  beam  about  8  feet  high  is  projected  onto  the  retro-reflective  screen  in 
order  to  detect  the  arrival  of  the  projectile  at  a  station.  When  the  projectile  enters  the  beam,  a 
trigger  is  generated  to  flash  the  spark-gap  light  source. 
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5.  RETRO-REFLECTIVE  SCREEN 

A  3M  No.  7610  screen,  made  up  to  tiny  glass  beads  about  0.002"  in  diameter,  is  used  as 
the  surface  on  which  a  shadow  of  the  projectile  is  cast.  The  shadowgraph  formed  allows  the 
camera  to  be  focused  on  a  known  object  (screen)  since  the  projectile  can  be  anywhere  in  the  12 
foot  width  of  the  range.  Since  the  screen  is  retro-reflective,  light  from  adjacent  stations  is  not 
reflected  to  the  wrong  camera.  Also,  the  screen  has  a  light  gain  over  that  of  a  Lambertian 
surface.  The  gain  is  dependent  on  the  viewing  angle  to  the  source.  A  gain  of  about  40  is 
obtained  in  the  Eglin  range  with  the  spark  separated  by  5"  from  the  camera  viewing  axis,  but  a 
gain  of  about  1600  is  possible  if  the  light  source  is  viewed  on  axis  with  the  camera.  The  light 
gain  of  the  screen  vs.  viewing  angle  is  shown  in  Figure  2. 
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Figure  2.  3M  Screen  Data 
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SECTION  III 

EARLY  PROTOTYPE  RESULTS 


Many  tests  of  the  vidicon  cameras  and  system  were  done  on  a  static  projectile  hanging  in  a 
simulated  station  at  Los  Alamos,  complete  with  screen,  spark-gap,  and  infrared  trigger.  Also, 
other  measurements  were  done  on  the  cameras  and  spark-gap,  and  some  image  processing  has 
been  done  for  image  enhancement  of  the  images  produced.  Some  of  the  results  are  reported  here. 

1 .  SYSTEM 

The  system  was  assembled  as  described  in  Section  II.  Many  non-shadowgraph  images 
were  obtained  using  overhead  room  lighting  and  timed  exposures.  The  cameras  were  adjusted 
for  optimal  image  quality  and  sensitivity.  Software  was  written  to  monitor  two  signals  (start  and 
stop)  for  image  acquisition,  and  then  to  store  the  images  in  memory.  Displayed  on  the  19" 
monitor  were  512  x  512  portions  of  the  2048  x  1024  images. 

2.  LIGHTING 

An  attempt  was  then  made  to  obtain  shadowgraph  images  using  the  spark-gap  as  the  light 
source,  with  discouraging  results.  No  image!  The  sensitivity  of  the  camera  was  measured, 
using  a  TV  Optoliner  operated  at  2856  “K  color  temperature  and  neutral  density  filters.  The 
measured  sensitivity  was  close  to  the  manufacturer's  data.  The  output  of  the  spark-gap  was  then 
measured  (Reference  2)  both  spectrally  and  for  total  energy  using  a  spectrometer  and  a 
wavelength  independent  joulemeter  for  pulsed  sources.  The  output  of  the  spark-gap  was 
determined  to  be  about  10  microjoules,  which  is  about  a  factor  of  500  less  than  that  determined 
in  the  initial  study  of  the  range  parameters.  Figure  3  shows  the  new  spectral  measurement  on  the 
spark-gap,  which  is  quite  different  than  the  initial  measurement  in  that  it  does  not  extend  nearly 
as  far  into  the  red.  This  probably  accounts  for  the  disparity  in  the  output  intensity. 


Figure  3.  Spark-Gap  Intensity  vs.  Wavelength 
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3.  RESOLUTION 


The  system  resolution  was  measured  (Reference  3)  using  a  TV  Optoliner  to  produce  an 
image  of  a  standard  TV  resolution  chart.  This  image  was  focused  directly  on  the  vidicon  tube  of 
the  camera,  to  get  resolution  independent  of  the  lens,  since  the  lens  contribution  could  be  readily 
determined.  The  results  of  the  measurements  show  a  resolution  of  about  0.82  Line  Pair/mm  at 
the  screen,  which  meets  the  requirements  of  the  system  determined  in  the  initial  study. 
Additionally,  another  test  showed  that  we  could  detect  and  locate  a  1  mm  diameter  pin  (used  for 
determining  spin  or  rotation  of  the  projectile)  which  was  another  requirement  of  the  range.  This 
indicates  that  the  initial  range  resolution  measurement  and  determination  of  the  electronic  image 
resolution  requirements  was  correct. 

4.  COST 

One  of  the  objectives  of  the  prototype  system  was  to  demonstrate  the  feasibility  of 
instrumenting  the  entire  range  at  Eglin  with  a  high  speed  electronic  imaging  system.  The  major 
benefits  of  doing  this  would  be  to  dramatically  reduce  the  time  from  a  test  to  when  results  would 
be  available.  This  could  literally  change  the  data  reduction  time  from  weeks  to  hours!  For  a  test 
which  requires  several  set-up  shots,  with  time  gains  on  each  shot,  the  benefits  are  compounded. 

The  cost  of  each  camera  was  about  $25,000  in  small  quantities.  It  was  estimated  that  for 
the  100  cameras  required  to  instrument  the  range,  the  cost  could  be  reduced  to  about  $10,000  per 
camera  (no  guarantee),  but  that  an  image  buffer  would  be  required  for  each  camera,  since  the 
image  processor  could  not  read  the  images  from  100  cameras  fast  enough  to  prevent  intolerable 
dark  current  contributions  to  image  noise.  So  a  cost  of  about  $15,000  per  camera  and  image 
buffer  was  estimated,  which  results  in  a  system  cost  that  exceeds  the  Eglin  budget  for  the 
Aeroballistics  Facility  conversion. 

5.  SUMMARY 

Results  of  the  first  measurements  on  the  prototype  system  were  very  disappointing.  Both 
from  a  performance  standpoint  and  from  a  potential  cost  standpoint.  The  idea  of  replacing,  or  at 
least  augmenting,  the  film  camera  system  is  a  good  one.  Higher  productivity,  shortening  the 
data  reduction  time  and  improving  repeatability  of  measurements  all  seem  to  be  possible  with  an 
electronic  shadowgraph  system.  Unfortunately,  cameras  with  cost  to  fit  the  Eglin  budget  and  to 
do  the  job  do  not  seem  imminent  within  the  next  few  years. 

Therefore,  at  this  stage,  it  seemed  that  unless  another  cost  effective  method  could  be  found 
to  make  Aeroballistics  measurements  to  the  same  accuracy  as  the  current  film  method,  the  project 
must  be  abandoned.  Also,  the  prototype  system  cannot  be  used  to  even  demonstrate  the 
feasibility  of  the  electronic  shadowgraph  method,  since  the  cameras  are  not  sensitive  enough  to 
be  used  with  the  existing  spark-gap  light  source. 


SECTION  IV 


PROJECT  RE-EVALUATION 


At  this  point  in  the  prototype  development  project,  it  was  decided  that  we  must  re-evaluate 
some  of  the  initial  decisions  and  project  direction.  Some  of  the  most  pressing  topics  to  be 
considered  were: 

1 .  COST 

The  Eglin  project  management  has  already  projected  that  a  cost  of  no  more  than  $10,000 
per  camera/image  buffer  is  the  budget  limit.  CCD  cameras  now  having  the  required  resolution, 
available  from  only  one  source,  and  image  memory  will  cost  at  least  twice  that  much.  Since 
costs  of  CCD  cameras  with  array  formats  larger  than  1024  x  1024  are  not  anticipated  to  be  driven 
lower  by  the  commercial  market,  we  do  not  think  that  the  cost  of  2048  x  1024  cameras  will  be 
reduced  appreciably  in  the  near  future. 

2.  CCD  vs.  VIDICON  CAMERAS 

There  are  several  reasons  to  prefer  CCD  cameras  over  vidicons  for  our  application.  First, 
the  geometric  distortion  in  CCD  cameras  is  practically  zero,  while  there  is  about  1  to  3  percent  in 
vidicons,  which  would  have  to  be  corrected.  Also,  CCDs  do  not  have  the  blooming  problem 
that  vidicons  have,  and  are  usually  smaller  and  more  rugged  than  vidicons. 

Additionally,  a  test  was  conducted  on  a  standard  CCD  camera  showing  that  it  did  have 
enough  sensitivity  to  be  used  with  our  spark-gap  light  source,  whereas  the  vidicon  did  not.  So, 
CCDs  are  to  be  chosen  as  the  preferred  technology  for  the  Eglin  cameras,  meaning  further 
vidicon  camera  experiments  would  not  be  productive,  since  noise  and  imaging  features  would  be 
different. 

3.  IMAGING  REQUIREMENTS 

The  initial  study  determined  the  resolution  capability  of  the  existing  film  system  and  then 
proposed  to  replace  it  with  a  comparable  electronic  imaging  system.  However,  those  stringent 
resolution  requirements  result  in  requiring  cameras  that  are  too  expensive.  Instead  of  just 
duplicating  the  film  system  capability,  we  looked  at  what  actually  is  needed  at  the  Eglin  range.  It 
was  determined  that  the  basics  of  what  was  needed,  was  to  very  accurately  determine  the 
absolute  position  of  the  projectile  and  its  orientation  (pitch,  yaw,  and  roll)  at  each  station.  The 
accuracy  requirements  of  those  measurements  was  next  determined. 

The  pitch  and  yaw  of  a  projectile  is  now  determined  geometrically  at  Eglin.  Since  most 
models  are  symmetrical  about  one  or  more  axis,  a  point  at  the  front  tip  and  at  the  tail  on  the  same 
axis  are  determined.  The  angle  of  the  line  connecting  these  points  is  the  angle  of  the  projectile. 
Since  the  images  are  actually  projections  of  an  object  on  a  shadowgraph,  three-dimensional 
geometry  is  required  to  determine  the  actual  angle  of  the  projectile  from  two  orthogonal 
shadowgraph  views  of  the  same  projectile.  The  accuracy  requirements  of  determining  the  angle 
of  the  projection  can  be  determined  from  the  range  requirements. 
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Figure  4  shows  a  model  projection  with  front  (F)  and  rear  (R)  points  selected,  and  a  line 
drawn  between  them.  A  grid  is  superimposed  to  represent  a  2048v  x  1024h  array  used  to 
acquire  the  image. 


Figure  4.  Model  Projection 

The  error  in  determining  the  actual  location  of  a  point  on  the  projectile  is  approximately 
+/-  0.5  pixel.  The  maximum  error  of  the  angle  of  the  line  connecting  the  two  points  F  and  R  is 
described  by  the  equation: 


Error  =  SIN"1- 
max 


(degrees) 


where  p  is  the  spatial  dimension  of  a  pixel  and  L  is  the  length  of  the  projectile.  For  the  Eglin 
range,  the  shadowgraph  screen  measures  2438.4  mm  (96")  high  x  1219.2  mm  (48")wide.  The 
1219.2  mm  dimension  at  the  screen  maps  onto  the  1024  pixel  dimension  of  the  array,  thus  the 
horizontal  spatial  dimension  of  each  pixel  is: 


1219.2  mm 


=  1.19 


The  vertical  spatial  pixel  dimension  is  the  same  as  the  horizontal. 

The  shortest  projectile  tested  and  measured  at  the  Eglin  range  is  about  80  mm  long  x  20 
mm  diameter.  Using  Equation  2  for  the  pixel  spatial  dimension  and  L  =  80  mm  in  Equation  1 . 
we  get  an  angular  error  of  +/-  0.853  degrees  for  the  shortest  projectile.  Keep  in  mind  that  this  is 
the  error  when  the  front  and  rear  points  are  located  with  no  error.  A  one  pixel  error  (operator 
error)  in  locating  the  front  and  rear  points  could  result  in  three  times  the  error  determined  above. 
However,  for  the  comparison  with  the  electronic  imaging  system  we  will  assume  that  there  is  no 
location  error.  The  angular  error  of  +/-  0.853  degrees  for  the  shortest  projectile  is  the  error  that 
the  electronic  imaging  system  must  achieve  to  match  the  2048v  x  1024h  array. 

4.  IMAGE  PROCESSING  REQUIREMENTS 

In  the  electronic  imaging  system,  noise  affects  the  accuracy  of  the  processing  algorithms, 
so  any  noise  minimization  will  be  helpful.  This  requires  Filtering  or  other  noise  reducing 
techniques  to  be  used,  and  places  a  great  emphasis  on  careful  selection  of  image  processing 
techniques  that  must  be  employed  to  provide  a  clean  image  while  not  affecting  measurement 
accuracy.  Also,  if  the  background  noise  can  be  eliminated  (by  image  processing),  then  the  entire 
image  can  be  processed  instead  of  having  to  Find  the  projectile.  This  can  greatly  simplify  any 
image  processing  which  automatically  finds  the  projectile  and  determines  its  location  and 
orientation. 


& 


SI 


SECTION  V 


REDUCING  RESOLUTION  REQUIREMENTS 


Two  primary  methods  have  been  investigated  to  reduce  the  resolution  requirements  to  allow 
using  a  lower  resolution  (and  thus  lower  cost)  camera  at  the  range.  These  methods  are:  finding 
centroids  (first  moment  of  inertia)  and  second  moments  of  inertia  of  areas.  An  example  of  using 
the  centroid  to  more  accurately  locate  an  object  is  demonstrated  with  the  reference  beads. 

One  of  the  tasks  of  any  imaging  system  at  the  Eglin  range  is  to  accurately  locate  reference 
beads  strung  on  Kevlar  lines  at  each  station,  as  shown  in  Figure  5,  in  order  to  locate  the 
projectile  shadowgraph  with  respect  to  the  beads,  and  thus  the  absolute  location  in  the  range. 
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Figure  5.  Shadowgraph  Station  Reference  Beads 

The  absolute  location  of  all  beads  in  the  range  are  accurately  known,  and  thus  provide  a 
reference  for  determining  the  absolute  location  of  the  projectile  shadowgraph,  since  both  the 
reference  bead.  Figure  6,  and  projectile  shadowgraph  are  contained  in  the  same  image. 
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All  Dimensons  in  inches 


Area  A  -  0.01472 
Area  S  -  0.023058 


Figure  6.  Eglin  Bead 
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The  projected  area  of  the  bead  on  the  screen  is  about  the  same  size  as  the  bead,  or  24.373 
mm2.  Since  the  bead  is  located  about  1/4"  from  the  screen,  there  is  very  little  magnification  of 
the  shadow  over  the  actual  bead  size.  The  bead  projected  area  covers  about  17  square  pixels  in  a 
2048  x  1024  array  (24.373  p2  =  17.2).  Locating  an  edge  of  the  bead  by  locating  a  point 
(without  further  image  processing)  would  still  have  a  +/-  0.5  pixel  error,  or  about  +/-  0.595  mm 
(0.0234")  error. 

However,  locating  the  centroid  of  the  bead  projection  on  the  shadowgraph  is  much  more 
accurate  than  the  geometric  method  described  above.  Each  pixel  of  the  bead  image  is  digitized  to 
12-bit  accuracy  (i.e.,  1  part  in  4096)  in  the  VICOM  system.  The  numerical  value  of  the  pixel 
(values  between  0  and  1)  represents  the  pixel  area  covered  by  the  image,  with  1  =  fully  covered. 
If  the  signal-to-noise  ratio  of  the  system  is  less  than  40  db,  the  noise  will  be  the  primary 
contributor  to  the  pixel  value  error,  and  the  digitizer  error  can  be  neglected.  A  signal-to-noise 
ratio  of  40  db  (typical  of  many  cameras)  is  a  voltage  ratio  of  100/1,  equivalent  to  a  pixel  area 
error  of  +/-  1  percent.  A  description  of  the  centroid  and  second  moment  of  inertia  algorithms, 
and  a  determination  of  overall  error  of  each  follows,  along  with  the  error  of  locating  the  bead 
from  the  example  above. 

1 .  CENTROID 

Figure  7  shows  a  general  case  for  determining  the  centroid  of  an  area  relative  to  a  set 
of  orthogonal  axes,  labeled  x  and  y: 


Y 


Figure  7.  Centroid  of  an  Area 


The  equations  (Reference  4)  for  the  centroid  of  the  area  A  are  given  by: 
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Where  AA  is  one  pixel  area,  fn  is  the  fraction  of  each  pixel  area,  AA,  covered  and  N  is  the 
number  of  pixels.  For  digital  systems,  the  integral  is  easily  approximated  by  doing  a 
summation. 

If  each  AA  is  accurate  to  +/-  1  percent  (due  to  noise  contribution),  then  the  sum  of  AAs  for 
17  pixels  (bead  example  with  a  2048  x  1024  array)  is  accurate  to  +/-  0.24  percent,  assuming 

noncorrelated  noise,  since  we  get  the  statistical  improvement  of  -7=  where  N  =  17.  The  location 

VN 

accuracy  of  0.24  percent  within  a  pixel  dimension  of  1.19  mm  yields  0.0029  mm  (0.0001") 
accuracy,  or  about  two  orders  of  magnitude  improvement  over  the  +/-  0.5  pixel  or  50  percent 
accuracy  computed  earlier  for  locating  a  point  Therefore,  we  see  that  the  centroid  method  of 
locating  an  object  is  more  accurate  than  just  locating  a  point  on  the  object,  provided  N>1. 

2 .  SECOND  MOMENT  OF  INERTIA  OF  AN  AREA 

Figure  8  shows  a  general  case  for  determining  the  second  moment  of  inertia  of  an 
area  with  respect  to  the  centroid  of  the  area.  It  turns  out  that  Finding  the  second  moment  with 
respect  to  the  centroid  rather  than  to  arbitrary  axes  is  simpler  because  the  arbitrary  axes  moment 
terms  do  not  have  to  be  carried  through  the  equations. 

AXIS  OF  MINIMUM  SECOND 
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The  equation  (Reference  4)  for  the  angle  of  the  axis  through  the  centroid  of  the  second 
moment  of  inertia  of  an  area  is: 


TAN  20  = 


h  -  V 


Where:  Iy  =  Jx'2dA 

Ix  =  jy'2dA 


(second  moment  with  respect  to  y  axis) 


(second  moment  with  respect  to  x  axis) 


Pc  =  jx’y'dA 


(product  of  inertia) 


and  0  =  Angle  of  the  axis  of  the  second  moment  of  inertia  relative  to  the  x'  axis, 

x'  =  Distance  from  the  centroid  to  dA  in  the  x  direction, 
y’  =  Distance  from  the  centroid  to  dA  in  the  y  direction. 

Note  that  the  result  yields  two  axes  90°  apart,  being  the  maximum  and  minimum  axes  of 
the  second  moment  of  inertia  of  the  area  A.  For  our  purposes,  we  are  only  interested  in  the  axis 
of  the  minimum  second  moment 

One  might  intuitively  surmise  that  finding  the  angle  of  the  second  moment  of  inertia  of  an 
area  is  more  accurate  than  the  geometric  solution  due  to  the  statistical  improvement  of  using  more 
data  points.  There  is  some  additional  help  due  to  the  fact  that  only  pixels  on  the  edge  of  the  area 
are  affected  by  noise.  Interior  pixels  are  known  to  be  completely  covered  by  the  image,  and  so 
their  values  can  be  corrected  (by  image  processing)  to  the  full  value  of  1.0. 

Rather  than  undertake  an  error  analysis  of  the  second  moment  of  inertia  for  random, 
uncorrelated  noise  on  edge  pixels,  a  simulation  was  used  to  prove  the  accuracy  of  the  technique. 

3 .  CENTROID/SECOND  MOMENT  OF  INERTIA  SIMULATION 

Using  the  equations  for  the  centroid  and  second  moment  of  inertia  of  an  area,  a 
general  case  problem  was  created  to  simulate  a  real  condition  to  demonstrate  the  accuracy  of  the 
method.  A  simulation  program  was  written  to  solve  for  the  second  moment  axis  and  use 
randomly  generated  "noise"  to  affect  the  edge  areas  of  the  projectile.  This  allowed  us  to 
determine  the  effects  of  random  errors  on  the  AA  values,  and  also  to  compute  the  standard 
deviation  of  our  calculations,  giving  us  a  feel  for  the  accuracy  we  might  expect  from  the  second 
moment  method. 


a.  Array  Size 


An  array  size  of  5 12v  x  256h  was  chosen  because  it  was  close  to  optimum  for 
maximum  size/cost  for  available  CCD  cameras,  at  the  time  of  the  simulation. 

b.  Projectile  Shape 

The  projectiles  generally  used  at  Eglin  are  bullets,  bombs,  and  shapes  that  fly 
well  without  tumbling  as  they  travel  downrange.  These  are  generally  oblong  shapes  with  the 
length  much  greater  than  the  diameter.  Also,  the  shapes  tend  to  be  symmetrical  about  the  long 
(horizontal)  axis.  To  make  matters  simple  for  the  drawings  and  simulation  area  calculations,  a 
bullet  shape  of  80  mm  long  by  20  mm  diameter  was  chosen,  because  this  is  the  the  smallest 
projectile  used  at  the  range.  This  casts  a  shadow  whose  area  is  easy  to  calculate.  For  the 
simulation,  the  projectile  was  drawn  at  an  angle  of  20  degrees  relative  to  a  horizontal  axis  (X). 

c .  Geometric  Model 

The  representative  model  used  is  shown  in  Figure  9  with  each  small  square 
representing  one  pixel  of  the  imaging  array.  All  units  of  measure  are  in  pixel  units. 


Figure  9.  Simulation  Model 

The  size  of  the  projectile  at  the  reflective  screen  was  then  mapped  onto  the  5 1 2V  x  256h 
array  using  a  spatial  resolution  factor  of  the  field  of  view  of  1219.2  mm  horizontal  +  256 
horizontal  pixels  =  4.76  mm  per  pixel  (same  for  the  vertical  dimension).  The  size  of  the  model 
in  pixel  units  is  16.8  pixels  long  by  4.2  pixels  high.  The  model  was  angled  at  20  degrees 
relative  to  the  horizontal. 

The  interpretation  of  the  shadowgraph  of  the  model  overlaying  the  array  is  that  it  simulates 
the  image  of  the  projectile  cast  on  a  CCD  array.  The  CCD  array  provides  an  output  to  the 
imaging  system  that  digitizes  to  12  bits  with  values  between  0.0  and  1.0.  Pixels  fully  covered 
by  the  shadow  will  have  a  value  of  1.0  in  the  system  with  no  noise  or  error  because  we  know 
they  are  interior  to  the  shadow.  These  are  the  pixels  shown  shaded  in  Figure  9.  The  pixels  that 
are  only  partially  covered  by  the  shadow  (edge  pixels)  will  output  a  value  between  0.0  and  1 .0. 
with  random  noise  added  to  the  ideal  value.  The  noise  is  a  random  value  between  0  and  +/-  the 
reciprocal  of  the  signal  to  noise  ratio  (S/N).  These  are  the  white  portions  of  pixels  shown  in 
Figure  9  that  are  inside  the  model  shadow  outline.  The  ideal  value  of  each  pixel  output  is 
proportional  to  the  pixel  area  subtended  by  the  shadow.  To  achieve  a  perfect  shadowgraph 
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l  model,  the  values  for  all  pixels  were  read  manually  from  Figure  9,  by  estimating  the  number  of 

■  subareas  on  a  Finer  gridded  graph  paper.  These  areas  were  used  in  the  simulator  as  the  baseline 

(data  for  pixel  areas.  The  areas  determined  manually  for  all  pixels  were  entered  in  the  simulator 
for  the  following  calculations,  with  random  noise  then  added  to  each  pixel  area. 

\ 

;  d.  Simulation  Results 


;  The  parameters  calculated  from  the  "perfect"  shadow  with  no  noise  added  are: 

i  Total  Area  =  59.171,  Cx  =  10.020,  Cy  =  3.562,  8  =  19.963  ° 

% 

K 

*  Table  1  shows  the  average  values  produced  by  the  simulator  with  noise  for  the  Area,  Cx,  Cy, 

*  angle  0,  standard  deviation  of  angle  0,  and  the  3a  value  for  angle  8 .  The  value  for  a  should  be 

I  used  as  the  probable  error  to  be  compared  with  the  current  probable  error  values  for  the  range 

p  using  the  film  cameras. 

i 

|r  The  S/N  ratio  is  calculated  using  the  root-mean-square  (RMS)  value  of  the  random  noise 

r!  generated,  which  is  a  continuous  function,  and  using  the  voltage  ratio  of  the  signal  and  noise. 

[■  The  signal  range  is  0  -  1  volt.  For  each  S/N  ratio,  the  simulator  was  run  25  times  to  provide  a 

I  good  statistical  sample.  The  standard  deviation  is  based  on  these  25  runs. 

* 

I 

l  TABLE  1.  SIMULATION  RESULTS  VS.  S/N 


S/N 

Cx 

Cy 

9 

a 

3a 

3.5 

10.031 

3.573 

19.082 

0.7931 

2.3793 

5.2 

10.007 

3.560 

19.995 

0.5614 

1.6842 

6.9 

10.055 

3.567 

19.782 

0.4198 

1.2594 

8.7 

10.021 

3.557 

19.879 

0.3596 

1.0788 

10.4 

9.991 

3.554 

20.054 

0.2658 

0.7974 

12.1 

10.011 

3.559 

19.899 

0.2357 

0.7071 

13.9 

10.016 

3.564 

19.984 

0.2389 

0.7167 

15.6 

10.018 

3.563 

20.001 

0.1920 

0.5760 

17.3 

10.017 

3.562 

19.984 

0.1225 

0.3675 

26.0 

10.016 

3.561 

19.974 

0.0992 

0.2976 

34.7 

10.020 

3.562 

19.950 

0.1021 

0.3063 

52.0 

10.016 

3.560 

19.968 

0.0617 

0.1851 

69.3 

10.020 

3.562 

19.964 

0.0399 

0.1197 

86.7 

10.020 

3.563 

19.967 

0.0199 

0.0597 

121.3 

10.021 

3.562 

19.964 

0.0192 

0.0576 

173.3 

10.020 

3.563 

19.969 

0.0172 

0.0516 

V.v.v.  V.  r. 
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Figure  10  shows  a  plot  of  the  3a  angle  0  error  vs  S/N  taken  from  Table  1.  It  also  shows 
equivalent  plots  for  5  and  10  degree  angles,  showing  insensitivity  to  model  orientation. 
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Figure  10.  3a  Angle  0  Error  vs.  S/N 
e.  Error  vs.  Array  Size 

Next,  in  order  to  look  at  the  effects  of  array  size  on  the  angular  error,  the 
simulator  was  modified  to  accommodate  various  array  sizes.  Simulations  were  run  for  2048  x 
1024,  1024  x  512,  and  512  x  256,  as  the  three  most  likely  array  sizes  to  be  used.  Figure  1 1 
shows  the  angular  error  (3a)  vs.  S/N  for  the  three  array  sizes.  These  plots  assumed  random 
noise  over  the  entire  projectile  image,  so  the  errors  in  general  are  higher  than  that  calculated  in 
paragraph  V.3.d  above. 


0  20  40  60  80  1  00  1  20 

S/N  RATIO 
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(1024  x  2048) 
(512 X 1024) 
(256  x  512) 


Figure  11.  Angle  Error  vs.  S/N  for  Various  Arrays 


P  ,  StN-I  -  SIN-1 


0.69  mm 
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The  P  error  of  0.494  degrees  is  based  on  the  assumption  that  an  operator  can  locate  the  end 
points  of  the  line  perfectly,  resulting  in  +/-  one  error  limit  of  resolution.  The  error  for  a  digital 
imaging  system  with  an  array  of  2048v  x  I024h  resolution  was  shown  to  be  +/-  0.852  degrees 
with  no  operator  error  in  locating  the  end  points  of  the  line.  The  error  would  be  three  times  that, 
or  +/-  2.56  degrees  if  there  were  a  +/-  1  pixel  error  in  finding  the  line  end  points.  Similarly  for 
the  film  system,  with  +/- 1  error  limit  in  mislocating  the  end  points  of  the  line,  a  1.48  degree 
error  would  result. 

By  comparison,  a  512v  x  256h  array  with  a  S/N  of  20,  which  is  realizable  with  the  Eglin 
reflective  screen  and  typical  CCD  cameras,  will  yield  an  angular  error  of  about  0.40  degrees.  It 
can  be  seen,  then,  that  using  the  second  moment  of  inertia  algorithm  with  a  5 12V  x  256h  CCD 
array  will  yield  angular  accuracies  of  about  what  the  present  film  system  can,  and  more  than 
twice  the  accuracy  of  a  2048 v  x  1024h  array  camera  using  the  geometric  method.  A  1024v  x 
512h  array  camera  then  will  yield  angular  determinations  better  than  the  film  system. 
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f.  Comparison  With  the  Film  Method 

The  Los  Alamos  National  Laboratory  report  #LA- 10072-MS  (Reference  5), 
contains  an  analysis  of  the  Eglin  Aeroballistics  Research  Facility.  This  analysis  shows  that  the 
resolution  capability  of  the  film  camera  system  used  at  the  Aeroballistics  Facility  has  a  limiting 
resolution  of  0.72  line  pair/millimeter  at  the  screen.  This  is  1.39  mm/line  pair  or  0.69  mm/line. 

This  means  that  an  object  0.69  mm  wide  can  just  be  discriminated,  or  0.69  mm  is  the  "error 
limit"  of  the  Eglin  film  system.  Using  that  number  for  the  error  dimension  for  locating  a  point 
on  each  end  of  a  line,  and  using  a  length  of  80  mm  for  the  shortest  projectile,  we  can  calculate 
the  error  in  finding  the  angle  of  a  line  on  film  by  using  Equation  1.  The  calculation  is  shown  in 
Equation  6. 
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SECTION  VI 


LOCATING  A  MODEL  IN  3-SPACE 


By  calculating  the  centroid  to  locate  an  object  and  the  second  moment  of  inertia  to  determine 
the  orientation  of  the  object,  the  location  and  orientation  can  be  determined  with  much  greater 
accuracy  than  by  using  geometric  methods  utilizing  only  one  or  two  points. 

This  section  details  the  algorithms,  equations,  and  terminology  for  determining  the 
orientation  of  a  projectile  model  and  the  location  of  the  center  of  mass,  Cg>  of  the  model  in  3- 
space  from  the  centroid  and  second  moment  of  inertia  data  from  two  orthogonal  electronically 
imaged  point  projections  taken  in  the  Eglin  Aeroballistics  Research  Facility  configuration. 

1 ,  GIVEN  INFORMATION 

a .  The  Eglin  Aeroballistics  Range  configuration  is  shown  in  Figure  12. 

b .  The  model  is  symmetric  about  its  principal  (horizontal)  axis. 

c .  The  model  is  imaged  on  orthogonal  Hall  (side  view)  and  Pit  (vertical  view) 
screens  by  point  flash  sources  at  F^  and  Fp,  with  coordinate  locations 

(Xfh’Yfh’  Zfh>  and  (Xfp’V  Zfp}  known- 
d  ,  The  electronic  imaging  cameras  have  detector  array  sizes  of  102^  x  5 1 2h  pixels. 

e  .  The  screen  size  (field  of  view)  to  be  imaged  is  243.84  cm  (96”)  high  x  121 .92  cm 

(48")  wide. 

f .  The  lower-right  calibration  bead,  B1Hi,  of  the  Hall  view  at  window  1  is  at 
range  coordinate  system  location  (0,0,0)  for  X,  Y,  Z. 

g  .  The  range  coordinate  system  has  the  X-axis  downrange,  not  necessarily  through 
B2H1.  the  Y-axis  crossrange,  and  the  Z-axis  vertical.  All  coordinates  are  given  in 
centimeters  and  are  referenced  to  this  coordinate  system. 

h  .  The  center  of  mass,  Cq,  as  well  as  the  length,  diameter,  and  centroid  of  the  side 
view  cross-sectional  area  are  known. 
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Figure  12.  Eglin  Aeroballistics  Range  Imaging  References 
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a .  Calculate  the  position  of  the  model  based  on  the  locations,  Ycg,  Z cg  of 

its  center  of  mass,  CG,  in  3-space  at  the  station  used. 

b  .  Calculate  the  orientation  of  the  principal  axis  of  the  model,  in  3-space,  by  the 

direction  cosines  (Reference  6)  COS  aw,  COS  Bw,  COS  yw,  of  the  model’s  principal 

M  M  M 

axis,  relative  to  the  three  range  coordinate  axes. 


3.  DEFINITIONS 
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The  following  designations  (Table  2)  will  be  used  as  reference  terms  in  the  descriptions  and 
equations  in  this  report  Designations  may  be  used  in  combination  Example:  C$H  refers  to  the 
(C)entroid  of  the  (S)hadowgraph  of  the  (H)all  view. 


TABLE  2.  REFERENCE  DESIGNATIONS 


R 

=  Spatial  resolution  of  image  (cm) 

S 

=  Shadowgraph  image  formed  on  the  Hall  or  Pit  screen 

c 

=  Centroid 

* 

B 

=  Bead 

0, 

5H’  0SP 

=  Angle  of  Hall  or  Pit  shadowgraph  relative  to  range  X  axis 

9, 

ISH’  9ISP 

=  Angle  of  Hall  or  Pit  shadowgraph  relative  to  image  h  axis 

9’ 

ISH’  6  ISP 

=  Angle  of  Hall  or  Pit  shadowgraph  relative  to  image  h'  axis 

H 

=  Hall  (view  on  wall  of  range) 

P 

=  Pit  (view  on  ceiling  of  range) 

I 

=  Image 

* 

F 

=  Flash  strobe  (This  is  a  point  light  source) 

A 

=  Point  calculated  where  the  axis  of  the  minimum  second  moment  of 
inertia  of  the  screen  projection  intersects  the  Z  axis  in  the  Hall 
view,  or  the  Y'  axis  in  the  pit  view 

★ 

Ml 

=  Model  length  (cm) 

* 

Md 

=  Model  diameter  (cm) 

* 

H:g 

=  Model  center  of  mass  relative  to  the  front  of  the  model  (cm) 

Me 

=  Model  centroid  relative  to  the  front  of  the  model  (cm) 

* 

E 

=  Coordinates  of  inside  edge  of  bead 

a,b,c 

=  Direction  numbers  of  a  line  in  3-space,  in  X,  Y,  Z  directions 

a.P.Y 

=  Direction  angles  from  the  X,  Y,  Z  axes  to  a  line  in  3-space 

* 

=  Known  quantity,  both  location  and  physical  dimensions 

4.  COORDINATE  SYSTEMS 


There  are  four  reference  coordinate  systems  used  in  the  calculations  for  determining  the 
location  of  the  model  from  two  orthogonal  digital  images. 

a.  Image  Coordinate  System  (h,v) 

Image  calculations  are  all  done  relative  to  the  image  coordinate  system  provided  by 
the  image  processing  software  (Reference  7)  running  on  an  IBM  PC/AT  computer.  The  image 
coordinate  system  is  shown  in  Figure  13. 


h  (orizontal) 


v  (ertical) 


Figure  13.  Image  Memory  Pixel  Address  References 

The  h  axis  is  positive  to  the  right  and  the  v  axis  is  positive  down.  All  dimensions  are  in 
pixels.  The  coordinate  location,  h  =  0,  v  =  0,  is  in  the  center  of  the  upper  left  pixel  of  the 
camera,  which  images  a  point  on  the  upper  left  of  the  shadowgraph  screen  when  viewed  from  the 
camera  position.  The  image  coordinate  range  is  5 1 2  pixels  in  both  the  h  and  v  directions, 
although  the  image  viewed  on  the  monitor  is  5 1 2^  x  480v.  This  represents  only  a  portion  of  the 
1024V  x  512^  image.  A  more  powerful  image  processing  system  is  required  to  process  the  full 
1024v  x512h  image.  The  image  coordinate  system  is  not  necessarily  parallel  to  any  other 
coordinate  system. 


b.  Bead  Coordinate  System  in  the  Image  (IT,  v') 

The  bead  coordinate  system  is  formed  in  the  image  by  placing  the  h'  axis  through  the 
centroids  of  beads  Bj  and  B2  (Cjb  1  and  Qb2  in  Figure  14).  The  vertical  axis  is  perpendicular  to 
IV  and  has  its  origin  at  Ci b  1  ■  They  are  designated  h',  v'  as  shown  in  Figure  14.  The  dimensions 
are  in  pixels.  This  bead  coordinate  system  is  not  necessarily  parallel  to  any  other  coordinate 
system. 
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Figure  14.  Bead  Coordinates  in  Image 

c.  Range  Coordinate  System  (X,  Y,  Z) 

The  range  coordinate  system  is  designated  X,  Y,  Z  and  is  shown  in  Figure  12.  The 
dimensions  are  in  cm.  Location  (0,0,0)  is  at  the  inside  edge  (not  the  centroid)  of  bead  Bim 
labeled  point  E  in  Figure  6.  The  X  axis  is  downrange  and  is  parallel  to  a  water  level  system  in 
the  range.  The  Y  axis  is  crossrange  and  is  also  parallel  to  the  water  level.  The  Z  axis  is  vertical 
and  perpendicular  to  the  X  and  Y  axes. 

d.  Station  Coordinate  System  on  the  Image  (X',  Y\  Z') 

A  two-dimensional  station  coordinate  system,  parallel  to  the  range  coordinate  system 
is  constructed  on  a  Hall  or  Pit  image  to  allow  transforming  image  coordinates  into  range 
coordinates.  Figure  16  shows  the  station  coordinate  system,  designated  X'|(,  Z'H,  in  the  Hall 
view,  and  X'p,  Z'p,  in  the  pit  view  (Q1  represents  either  Z’H  or  Z’p).  The  origin  of  each  is  at  th 
centroid  of  bead  Bp  These  station  coordinate  axes  are  simply  offset  from  the  range  axes  by 
constants,  which  are  known,  because  the  locations  of  the  beads  in  the  ranee  are  known. 


5.  IMAGE  CALCULATIONS 


The  image  calculations  will  produce  the  centroid,  Cjs,  of  the  shadowgraph  of  the  model.  A 
subscript  "I"  is  used  to  denote  references  to  image  calculations,  relative  to  either  the  h,v  or  h',v' 
axes.  A  '(prime)  is  used  to  distinguish  between  the  h,v  axes  and  the  h',v’  axes.  The  centroid 
units  are  in  pixels  in  the  h,v  directions.  Also  produced  is  the  minimum  second  moment  of 

inertia,  9  [5,  through  the  centroid  of  the  shadowgraph,  created  by  illuminating  the  model  from  a 
point  light  source,  as  shown  in  Figure  12.  The  angle,  0IS,  of  the  second  moment  axis  relative  to 
the  h  axis  is  positive  in  the  positive  v  direction,  so  that  0°  lies  parallel  to  the  h  axis  and  +90°  is 
down,  parallel  to  the  v  axis  (Figure  13).  Then  the  centroids  and  second  moments  are  translated 
to  the  bead  reference  axes,  preparatory  for  transformation  to  the  range  coordinate  system. 

a.  Hall  Shadowgraph  Image  Calculations 

( 1 )  Calculate  the  centroid,  Cjbih  (hicB  1H»  vICB1h)>  °f  bead  Bjpj  and  the  centroid, 
CIB2H  (hICB2H>  vICB2h)>  of  bead  B2H- 

(2 )  Calculate  the  centroid,  (bicsH>  vICSh)>  °f  the  shadowgraph  of  the  model 
on  the  Hall  screen. 

(3 )  Calculate  the  angle,  G^h  (degrees),  of  the  minimum  second  moment  of  inertia 
of  the  shadowgraph  area,  through  the  centroid  CISH,  relative  to  the  hH  axis. 

b.  Pit  Shadowgraph  Image  Calculations 

( 1 )  Calculate  the  centroid,  Cjbip  (hjQjjp,  vjCB1p),  of  bead  B1P  and  the  centroid, 
CIB2P  (hiCB2P>vlCB2p)>  bead  B2p. 

(2 )  Calculate  the  centroid,  Cjsp  (hjcsP>  VICSP)>  °f  the  shadowgraph  of  the  model 
on  the  Pit  screen. 

(3 )  Calculate  the  angle,  9jsp  (degrees),  of  the  minimum  second  moment  of  inertia 
of  the  shadowgraph  area,  through  the  centroid  CiSP,  relative  to  the  hP  axis. 

c.  Transform  Image  Coordinates  to  the  Image  Bead  Coordinate  System 

Figure  15  shows  the  bead  coordinate  axes,  h',  v',  relative  to  the  image  axes,  h,v. 

The  transformation  can  be  derived  once  to  apply  to  both  the  Hall  and  Pit  images,  with  all 
dimensions  in  pixels.  To  obtain  equations  for  the  Hall,  add  the  subscript  H;  for  the  Pit,  add  the 
subscript  p.  The  points  and  angles  to  be  transformed,  which  are  calculated  by  the  centroid  and 
second  moment  of  inertia  methods,  are  indicated  as  follows  to  simplify  the  notation  for  the 
derivation  of  all  coordinate  transformations: 

Point  1  =  Centroid  of  Bead  Bj 

Point  2  =  Centroid  of  Bead  B2 

Point  3  =  Centroid  of  projectile  shadowgraph 

0j  =  Angle  of  minimum  second  moment  of  inertia  of  shadowgraph  relative  to  h 

9'j  =  Angle  of  minimum  second  moment  of  inertia  of  shadowgraph  relative  to  h' 


y*jr*ji  v  mj  w  vw  v  v  viv-v.yj#, 


Figure  15.  Image  to  Bead  Axis  Transformation 

Locations  of  the  points  in  the  image,  relative  to  the  h,v  axes,  are  designated: 

Point  1  =  hi,  v j  =  Cjbi  =  hrcB1,  vICBi  (7) 

Point  2  =  h2,  v2  =  CiB2  =  hlCB2>  V1CB2 

Point  3  =  h3,  v3  =  CIS  =  hics,  vICS  (9) 

(1)  Locate  the  Bead  Coordinate  System  in  the  Image 

First,  the  h'  axis  is  drawn  through  Point  1  and  Point  2  with  the  origin  located  at 
Point  1,  as  shown  in  Figure  15.  Then  the  v'  axis  is  constructed  perpendicular  to  the  h'  axis, 
with  its  origin  at  Point  1.  Point  1  will  now  become  the  origin  (0,0)'  for  the  h',v'  axes.  The 
dimensions  of  both  the  h,v  and  h',v'  axes  are  in  pixels. 

(2)  Compute  the  Angle 

The  angle  O  represents  the  rotation  between  the  h,v  and  h',v'  axes  and  is 
formed  by  the  h'  axis  and  a  line  through  Point  1,  parallel  to  the  h  axis.  The  triangle  formed  with 
sides  r,  e,  t  is  used  to  find  the  angle  O  which  is  positive  in  the  direction  shown.  Sides  r  and  e  are 
calculated  from  the  locations  of  Points  1  and  2,  which  are: 

r  =  vj  -  v2  (10) 

e  =  hi  -  h2  (11) 
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Then  the  angle,  C>,  is  given  by: 


O  =  TAN'1  -  =  TAN'1 -  (12) 

'  h,  -  h, 

(3)  Transform  the  Image  Points 

Because  the  centroids  for  Points  1,  2,  3  have  been  calculated  in  paragraph 
VI.5.a  -  VI.5.b,  equations  for  the  following  lines  (magnitude)  and  angles  in  Figure  15  are: 


j  =  h!  -h3 

03) 

k  =  vj  -  v3 

(14) 

g  =  ^J2  +  k5 

05) 

• 

.  ^r!,e5 

(16) 

0)  =  TAN-1  £ 

07) 

6  =  co  +  <P  =  TAN'1  £  +  TAN'1  T- 

08) 

So  the  locations  of  Points  1,  2,  3  and  angle  8'i  in  the  bead  coordinate  system  (h',v')  are: 

Point  1: 

h'j,  v'i  =(0,0)  (by  definition) 

09) 

Point  1: 

h'2,  V2  =  (t,  0) 

(20) 

Point  3: 

h'3,  v'3  =  (g  SIN  5,  g  COS  5) 

(21) 

0'l  =  0[  - 

(22) 

d.  Transform  the 

Bead  Coordinate  System  to  Station 

Coordinate  System 

The  equations  for  the  station  coordinate  system  can  be  derived  once,  using  station 
axes  X’,  Q’,  where  Q’  is  a  dummy  variable  to  be  replaced  by  either  Z'  for  the  Hall  station  or  Y' 
for  the  Pit  station  in  the  final  equations.  Likewise,  Q  is  a  dummy  variable  used  with  X  when 
referring  to  a  range  plane,  which  is  either  the  X,  Z  plane  or  the  X,  Y  plane.  Axis  Q  will  be 
replaced  by  Z  or  Y  when  writing  final  equations  for  the  Hall  or  Pit. 
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(1)  Transform  Bead  Coordinates  to  Station  Coordinates 

Figure  16  shows  the  reference  axes  X',  Q'  drawn  in  the  same  view  as  the  h',v' 
axes.  The  station  axis  X'  has  its  origin  at  the  centroid  of  Bead  Bj,  and  is  parallel  to  the  range  X 

axis.  The  station  Q'  axis  is  parallel  to  either  the  range  Y  or  Z  axis,  with  its  origin  at  the  same 
location  as  X'.  Point  locations  are: 

Point  1  =  Centroid  of  Bead  B 1 
Point  2  =  Centroid  of  Bead  B2 
Point  3  =  Centroid  of  projectile  shadowgraph 

Point  4  =  Intersection  of  second  moment  of  inertia  of  shadowgraph  with  Q'  axis 
(This  is  called  Point  A$  in  later  sections) 


h' 


Figure  16.  Bead  to  Station  Axis  Transformation 

The  origin  of  axes  X',  Q'  lies  on  the  origin  of  h',  v'  in  the  image.  The  dimensions  of  the 
X',  Q'  axes  are  in  pixels,  and  can  be  converted  easily  to  centimeters  by  knowing  the  spatial 
resolution  of  the  image. 

(2)  Compute  the  Angle  y 


The  angle  V|/  is  the  rotation  angle  between  the  h',v'  and  X',  Q'  axes.  The  range 
coordinates  for  the  edges  of  the  beads  are  known  as: 

Edge  of  Bead  Bj  =  XEbi>  YEbi>  Zebi  (23) 

Edge  of  Bead  B2  =  XEb2,  YEb2- Zeb2  (24) 

and  the  distance  from  the  edge  of  the  bead  to  its  centroid  is  shown  in  Figure  6  and  is  measured  as 

dB  =  0.2728  cm  (25) 
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Figure  17  shows  the  arrangement  of  the  beads  in  both  the  Hall  and  Pit  views  and  the 
relationship  of  the  centroids  of  the  beads  to  the  inside  edges,  which  are  known  in  the  range 
coordinate  system. 
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Figure  17.  Bead  Locations  in  Hall  and  Pit  Views 


Then  the  approximate  range  coordinates  of  the  centroid  of  the  beads  are: 

Cbi  =  (Xebi  '  d-B)>  ^EBl’  ZEB1  =  Yj,  Zj  (26) 

Cb2  =  (XEB2  +  de).  Yeb2>  ZEB2  =  X2>  Y2,  Z2  (27) 

Equations  26  and  27  are  not  exact  if  there  is  appreciable  sag  in  the  kevlar  line  connecting  the  two 
beads.  But  until  further  investigation,  these  equations  will  be  used  for  now. 

Since  the  X,  Q  and  X',  Q'  axes  are  parallel,  the  X,  Q  range  coordinates  for  Points  1  and  2 
can  be  used  to  compute  the  rotation  angle,  \| f,  as  shown  in  Figure  18. 


Figure  18.  Rotation  Angle  Calculation 
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The  location  of  Point  2  in  station  coordinates,  X',  Q',  is: 

X’^cm)  =  X2-Xj  (cm) 

Q'2(cm)  =  Q2  -  Ql  (cm) 


So  the  rotation  angle  is: 


y  =  TAN-1  =  TAN’1  ^ 


(3)  Transform  Image  Points  to  Station  Coordinate  System  (X',  Q') 

By  definition,  the  location  of  Point  l,  with  respect  to  to  X',  Q'  is: 

X'j.Q'p  (0,0)  (31) 

Knowing  the  angle  \|/,  we  can  locate  Point  2,  in  pixels,  with  respect  to  X’,  Q’  as: 

X'2  =  h’2COS\j/  (32) 

Q’2  =  h’2  SIN  \|/  (33) 

The  angle  5  is  shown  in  Equation  18.  From  Figure  16,  angle  a  is: 

a  =  (8  +  \|/)  (34) 

Then  the  location  of  Point  3  can  be  calculated  from  the  triangle  formed  by  sides  g,  X  3, 

Q'3,  where  g  has  been  calculated  in  Equation  15.  Point  3  is: 

X’3  =  gSINa  (35) 

Q'3  =  g  COS  a  (36) 

The  angle  of  the  second  moment  of  inertia  of  the  shadowgraph  is  simply: 

0’  =  0'j  -  v  (37) 

(4)  Locate  Point  4  (A§)  on  the  Q'  Axis 

The  location  of  Point  4  is  the  intersection  of  the  second  moment  of  inertia  of  the 
shadowgraph  axis  with  the  Q'  axis,  so  the  X'4  dimension  is  zero,  and  can  be  calculated: 


X'4 

=  0  (by  definition) 

(38) 

Q’4 

=  Q’3  *  n 

(39) 

n 

=  X'3  TAN  0’ 

(40) 

So: 

Q’4 

=  Q'3  -  X’3  TAN  0' 

(41) 

(5)  Transform  Hall  Bead  Coordinates  to  Station  Coordinates 

The  equations  derived  in  paragraphs  VI.5.d.(l)  -  VI.5.d.(3)  can  be  written  for 
the  Hall  image  Points  1,  2,  3,  4,  by  substituting  Z'h  for  Q’  terms: 


Point  1: 

x’h 

z'h 

' 

0 

0 

(42) 

(43) 

Point  2: 

X2H 

= 

h'2H  COS  V)/H 

(44) 

Z  2H 

= 

h'2H  SIN  yH 

(45) 

Point  3: 

X 

X 

= 

gH  SIN  aH 

(46) 

Z  3H 

= 

gH  COS  aH 

(47) 

Point  4: 

X 

■u" 

X 

= 

0 

(48) 

Z4H 

= 

Z3H  *  X3H  TAN  0'H 

(49) 

Line 

TH 

= 

V1H  -  V2H 

(50) 

Line 

eH 

= 

hlH  -  h2H 

(51) 

Line 

JH 

= 

h1H-h3H 

(52) 

Line 

kH 

= 

V1H  ‘  V3H 

(53) 

Line 

§H 

= 

+  kH 

(54) 

Line 

*H 

= 

h'2H  =  a/4  +  rH 

(55) 

Angle 

9’h 

= 

X 

> 

1 

X 

CD 

(56) 

Angle 

= 

(  5h  +  Vh) 

(57) 

Angle 

5h 

= 

“H  +  ^H 

(58) 

Angle 

coH 

= 

TAN’1  j*- 
kH 

(59) 

Angle 

oH 

= 

TAN'1  — 
eH 

(60) 

Angle  \|fH 


_  TAN-i  ZCB2H  ~  ZCB1H 
XCB2H  '  XCB1H 


(61) 


(6)  Transform  Pit  Bead  Coordinates  to  Station  Coordinates 


The  equations  derived  in  paragraphs  VI.5.d.(l)  -  VI.5.d.(3)  can  be  rewritten 
for  the  Pit  image  points  1,  2,  3,  4,  by  substituting  Y'p  for  Q'  terms: 


Point  1: 

X'p 

= 

0 

Y’p 

= 

0 

Point  2: 

X'2P 

= 

hSp  COS  tj/p 

Y’2p 

- 

h'2p  SIN  \\fp 

Point  3: 

X’3P 

= 

gp  SIN  oP 

Y  3p 

= 

gp  COS  Op 

Point  4: 

X'4P 

- 

0 

Y  4p 

= 

Y’3p  -  X'3p  TAN  0'p 

Line 

rP 

= 

V1P  ‘  V2P 

Line 

ep 

= 

hip-  h2P 

Line 

jp 

= 

hlP '  h3P 

Line 

kP 

= 

Vip-  v3P 

Line 

SP 

= 

^jiP  +  kp 

Line 

*P 

= 

1  2  2 
h  2P  =  +  rp 

Angle 

0'p 

= 

9  ip  -  Vp 

Angle 

Op 

= 

( 8p  +  \|tp) 

Angle 

5P 

= 

Up  +  Op 

Angle 

CDp 

= 

TAN'1 

kP 

Angle 

Op 

= 

TAN'1  T-£- 

Cp 

Angle 

Vp 

= 

tan-i 

aCB2P  ‘  *CB 

(62 

(63 


(64 

(65 


(66 

(67 


(68 

(69 


(70 

(71 

(72 

(73 


(74 

(75 


(76 

(77 

(78 


(79 

(80 

(81 
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e.  Transform  the  Station  Coordinates  to  Range  Coordinates 

This  translation  does  not  involve  any  geometry,  but  simply  requires  multiplying  X',  Y\  Z' 
locations  (in  pixels)  by  the  spatial  resolution  of  the  Hall  or  Pit  view,  and  then  adding  the  range 
location  offset  to  the  origin  of  the  X',  Y'  or  X',  Z'  station  views.  The  range  location  offset  to 
the  origin  of  the  Hall  or  Pit  station  is  the  location  of  the  Bead  B  j  centroid  in  each  case. 

( 1 )  Calculate  the  Hall  and  Pit  Station  Spatial  Resolutions 

The  spatial  resolution,  Rj,  of  any  image  is  defined  as: 


Rt  = 


distance  (cm) 
#  of  pixels 


The  distance  between  Bead  Bj  and  Bead  B2,  in  pixels,  has  been  calculated  in  each  view  as  the 
value  h'2  =  t  (pixels)  in  Equation  20.  The  range  distance  between  beads  along  the  h’  axis  for  the 
Hall  and  Pit  views  is  shown  in  Figure  18  and  can  be  calculated  in  X,  Y,  Z  units  (cm).  The  Hall 
and  Pit  spatial  resolutions  are  calculated  as  follows: 


Hall  Spatial  Resolution  in  cm/pixel  is: 


h'2H  (cm) 
h'2H  (Pixels) 


The  value  h'2H(PiXels)  ^as  b een  calculated  as  Equation  55.  The  value  for  h'2H(cm)  is  given  by: 

h2H(cm)  =  V(X'2H(cm))2  +  (Z'2H(cm))2  (cm)  (84; 

where  X'2H(cm)  is  from  Equation  28  and  Z'2H(Cm)  *s  fr°m  Equation  29,  substituting  Z'  for  Q' 
and  Z  for  Q.  Values  for  X^  X2  and  Z\,  Z2  are  known. 

The  pit  Spatial  Resolution  in  cm/pixel  is: 

r  -  h2P  (cm)  cm 

p  h’2p  (pixels)  pixel  ^ 

The  value  h'2p(pixels)  has  been  calculated  from  Equation  75.  The  value  for  h'2p(cm)  is  given  by: 

h’2H(cm)  =  V(X'  2H(cm))2  +  (z2H(cm))2  (cm)  (86: 

where  X'2p(Cm)  is  from  Equation  28  and  Y'2p(Crn^  is  from  Equation  29,  substituting  Y’  for  Q' 
and  Y  for  Q.  Values  for  Xj,  X2  and  Yj,  Y2  are  Known. 


(2)  Hall  Station  to  Range  Coordinate  Transformation 


The  final  equations  for  the  Hall  station,  in  terms  of  range  coordinates,  will  be 
obtained  by  multiplying  X'H,  Z'H  coordinates  by  and  adding  the  range  offset  tc  the  centroid 
of  Bead  Bj.  The  Hall  equations  in  paragraph  VI.5.d.(5)  for  lines,  with  dimensions  in  pixels, 
and  angles  are  still  valid  and  do  not  need  to  be  converted.  All  equations  will  also  show  the 
proper  subscripts  and  designations  for  the  points  and  angles. 


CB1H 

XCB1H’YCB1H>  ZcBlt^  (XEB1H  '  dB)<  YEB1H>  ZEB1H 

(87) 

CB2H 

=  XCB2H>  yCB2H»  ZCB2H=  (XEB2H  +  dB)«YEB2H>zEB2H 

(88) 

X  2H(cm) 

=  (XCB2H  -  XCB1H) 

(89) 

Z2H(cm) 

=  (ZCB2H  ‘  ZCB1H) 

(90) 

rH 

=  VCB1H  '  VCB2H 

(91) 

eH 

=  hCBlH  *  hCB2H 

(92) 

JH 

=  hCBlH  '  hCSH 

(93) 

kH 

=  VCB1H  ‘  VCSH 

(94) 

SH 

=  +  kH 

(95) 

lH 

=  th  +  eH 

(96) 

0>H 

JH 

=  TAN-1  ~ 

(97) 

oH 

tr 

=  TAN'1-— 
eH 

(98) 

0'lH 

i 

E 

<x> 

ii 

(99) 

5H 

=  C0h  +  Oh 

(100) 

Vh 

Z' 

2H(cm) 

y  * 

2H(cm) 

(101) 

=  SH  +  v)/H 

(102) 

0’h 

=  0'lH  +  Vh 

(103) 
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(103) 


But  0$h  =  0’n  since  the  X  and  X'  axes  are  parallel,  so: 


0sh 

=  0IH  ■  -  VH 

(104) 

CSH: 

XCSH 

=  RH  (gH  SIN  OH)  +  XCB1H 

(105) 

ycsh 

=  ycbih 

(106) 

ZcSH 

=  rh  (8h  COS  ch)  +  Zcbih 

007) 

Ash; 

xash 

=  XCBih 

(108) 

yash 

=  YCB1H 

(109) 

ZASH  =  (ZcSH  '  XCSH  TAN  0Sh)  +  ZCB1H  0 10) 

(3)  Pit  Station  to  Range  Coordinate  Transformation 

The  final  equations  for  the  Pit  station,  in  terms  of  range  coordinates,  will  be 
obtained  by  multiplying  X'p,  Y'p  coordinates  by  Rp  and  by  adding  the  range  offset  to  the 
centroid  of  Bead  Bp  The  Pit  equations  in  paragraph  VI.5.d.(6)  for  lines,  with  dimensions  in 
pixels,  and  angles  are  still  valid  and  do  not  need  to  be  converted.  All  equations  will  also  show 


the  proper  subscripts  and  designations  for  the  points  and  angles. 

CBip  =  XCB1P,YCB1P,  (XEB1P  -  dB),  YEBjp,  ZEB]P  (111) 

CB2P  =  XCB2P,  Ycb2P,  ZCB2p=  (XEB2p  +  dB),Y£B2P’ZEB2P  (H2) 

X'2P(cm)  =  (XCB2P  -  XCB1P)  (113) 

Z2P(cm)  =  (ZCB2P~  ZCBIP)  (114) 

rP  =  VCB1P  *  VCB2P  (1 15) 

eP  =  hCBlP'hCB2P  (116) 

JP  =  hCBlP '  hCSP  (1 17) 

kP  =  VCB1P  ’  VCSP  (1 18) 


gp 

=  A^Jp  +  kp 

(119) 

lP 

-  -  4 

(120) 

COp 

jp 

=  TAN-1  — 

(121) 

kP 
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V  w.  v  v  v-v.v» 


Op 

=  TAN’1— 
eP 

(122) 

0'lP 

=  0  jp  -  Op 

(123) 

SP 

=  cop  +  Op 

(124) 

Vp 

Z' 

T  AH-1  2P(cm) 

^  2P(cm) 

(125) 

aP 

—  8p  4-  \|Tp 

(126) 

e'p 

=  0'ip  +  Vp 

(127) 

But  0sp  =  0'p  since  the  X  and  X'  axes  are  parallel,  so: 

0SP 

=  0jp  -  Op  -  \|/p 

(128) 

CSP:  XCSP 

=  Rp  (gp  SIN  Op)  +  X^Bjp 

(129) 

YCSP 

=  RP  (gp  COS  aP)  +  YCB1p 

(130) 

ZCSP 

=  ZCBIP 

(131) 

ASP;  XASp 

=  XCB1P 

(132) 

YASP 

ZASP 

=  (YCSP  '  xcsp  tan  0sp)  +  YCbip 
=  ZCBIP 

(133) 

(134) 

6 .  LOCATE  THE  MODEL 

Using  the  information  from  the  shadowgraphs,  the  location  and  orientation  of  the  model,  in 
3-space,  is  determined.  This  is  accomplished  by  constructing  a  Hall  planen  and  a  Pit  planep 

using  the  shadowgraph  data  from  paragraph  VI.5.  These  planes  are  shown  in  Figure  12.  The 
intersection  of  the  two  planes  is  a  line,  which  is  the  principal  axis  of  the  model. 


a.  Determine  the  Planes 
(1)  Plane  Notation 

Referring  to  Figure  12,  Plane^  is  formed  by  three  points  Fh,  Csh-  Ash-  Planep  is  formed 
by  points  Fp,  C$p,  ASp.  The  points  have  been  defined  and  calculated  in  earlier  sections. 
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For  simplicity  sake,  all  points  and  angles  used  in  model  locating  calculations  now  will  be 
referred  to  as  follows: 

Fh  =  Hall  Flash 

Ch  =  Hall  shadowgraph  centroid 

Ah  =  Hall  shadowgraph  second  moment  axis  intersection  with  Z'  axis 

FP  =  Pit  Flash 

Cp  =  Pit  shadowgraph  centroid 

Ap  =  Pit  shadowgraph  second  moment  axis  intersection  with  Y'  axis 
where  the  locations  of  the  points  have  been  determined  as: 


Fh 

=  XFH>  YFH>  ZFH 

(135) 

CH 

=  XCSH>  YCSH>  ZCSH 

(136) 

ah 

=  XASH>  YASH>  ZASH 

(137) 

FP 

=  XFP>  Ypp,  Zpp 

(138) 

CP 

=  XCSP>  YCSP>  Zcsp 

(139) 

AP 

=  XASP»  YASP>  ZASP 

(140) 

The  angles  calculated  using  the  second  moment  of  inertia  method  are  designated: 

eH 

0P 

=  0SH 
=  0SP 

(141) 

(142) 

(2)  Derive  the  Plane  Equations  -  Planen  and  Planep 

The  two  planes  intersect  along  the  principal  axis  of  the  model.  The  equation  of 
a  plane  passing  through  three  points  (Reference  6)  is  given  by  the  form: 

JX 

+  KY  +  LZ  +  M  =  0 

(143) 

where  J,  K,  L,  and  M  are  constants  and  X,  Y,  Z  are  point  coordinates. 

Since  each  point  on  the  plane  must  satisfy  Equation  143,  then  three  equations  from  the 
three  points  may  be  written: 

JX! 

+  KY!  +  LZj  +  M  =  0 

(144) 

JX2 

+  KY2  +  LZ2  +  M  =  0 

(145) 

JX3 

+  ky3  +  lz3  +  M  =  0 

(146) 
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Solving  Equations  144,  145,  146  for  the  coefficients  J,  K,  L,  M  yields  the  determinants: 


YI 

Zi  i 

X1 

Zl 

1 

J  = 

Y2 

Z2  1 

(147) 

K  =  - 

X2 

Z2 

1 

Y3 

Z3  1 

X3 

z3 

1 

X1 

Y1  1 

X1 

Y1 

Zi 

L  = 

X2 

Y2  1 

(149) 

M  =  - 

X2 

Y2 

Z2 

X3 

Y3  1 

X3 

Y3 

z3 

(148) 


(150) 


If  the  three  points  defining  PlaneH  are  FH,  CH,  AH,  and  FP,  CP,  AP  define  PlaneP,  then  we 
will  let  Equation  144  be  the  equation  for  point  F,  Equation  145  represent  C,  and  Equation  146 
represent  A.  Solving  the  three  equations  144,  145,  146  for  the  coefficients  JH,  KH,  LH,  MH,  we 
obtain  the  equation  for  PlaneH,  which  is: 


JpjX  +  K^jY  +  Lj^Z  +  Mpj  =  0 

where  the  coefficients  J^,  KH,  LH,  MH  are  defined  by  the  determinants: 


yfh 

zfh 

1 

XFH 

ZFH 

l 

ych 

ZCH 

1 

(152) 

X 

ii 

i 

XCH 

ZCH 

l 

yah 

Zah 

1 

XAH 

zah 

l 

XFH 

YFH 

1 

XFH 

YFH  Zfh 

XCH 

ych 

1 

(154) 

mh=  - 

XCH 

ych  Zch 

XAH 

Yah 

1 

XAH 

Yah  Zah 

(151) 


(153) 


(155) 


Likewise,  solving  equations  144,  145,  146  for  JP,KP,LP,MP,  the  equation  for  Planep  is: 

JPX  +  KPY  +  LPZ  +  MP  =  0  (156) 

where  the  coefficients  JP,  KP,  LP,  MP  are  defined  by  the  determinants: 


Yfp  ZFP  1 

XFP  Zpp  1 

JP  = 

YCP  Zcp  1 

(157)  KP  =  - 

XCP  Zcp  1 

yap  zap  l 

x ap  ZAp  l 

LP  = 


XFP  YFP  1 
XCP  YCP  1 
XAP  YAP  1 


XFP  YFP  ZFP 
XCP  YCP  Zcp 

x ap  Yap  Zap 


(159)  MP  =  - 


(160) 


(3)  Compute  the  Intersection  of  the  Two  Planes,  Plane^,  Ptanep. 

The  intersection  of  the  two  planes  is  a  line  along  the  principal  axis  of  the  model. 
The  line  is  described  by  direction  numbers  aM,  bM,  cM,  or  by  the  direction  cosines  COS  CtM, 
COS  (3m,  COS  yM. 

Using  Equations  151  and  156,  the  direction  cosines  of  the  line  of  intersection  of  the  two 
planes  may  be  determined.  From  the  intersection  of  two  lines  lying  in  the  two  planes,  we  can 
determine  a  point  on  the  line  of  intersection  of  the  two  planes.  One  point  and  the  direction 
cosines  of  a  line  are  sufficient  to  fully  describe  the  line.  The  point  corresponds  to  the  centroid  of 
the  side  view  cross-sectional  area  of  the  model,  because  that  is  what  is  determined  from  the  Hall 
and  Pit  images.  The  line  corresponds  to  the  principal  axis  of  the  model,  since  that  is  determined 
also  in  the  images.  The  coefficients  of  X,Y,Z  in  the  equation  for  a  plane  are  proportional  to  the 
direction  cosines  of  a  line  perpendicular  to  the  plane.  Because  the  line  of  intersection  of  two 
planes  must  lie  in  both  planes,  then  it  must  also  be  perpendicular  to  the  perpendicular  lines  of 
both  planes  described  by  their  direction  cosines.  Because  the  direction  cosines  of  a  line  are  also 
proportional  to  the  direction  numbers  of  the  line,  direction  numbers  may  be  used  to  define  a  line. 

The  requirement  for  perpendicularity  for  two  lines  (Reference  6)  is: 

a  a1  +  b  b1  r  c  c'  =  0  (161) 

where  a,b,c  are  the  direction  numbers  for  one  line,  and  a',b',c'  are  the  direction  numbers  for  the 
second  line. 

The  coefficients  of  X,  Y,  Z  from  Equations  151  and  156  are  J,  K,  L,  and  they  are 
proportional  to  the  direction  numbers  for  the  line  perpendicular  to  the  plane  that  describes  the 
plane.  Substituting  the  coefficients  for  Equations  151  and  156  for  the  direction  numbers  a',b',c' 
in  the  previous  equation,  two  equations  result: 

Jpj  a  +  Kpj  b  +  Lh  c  =  0  (162) 

Jpa  +  Kp  b  +  Lpc  =  0  (163) 

Solving  Equations  162  and  163  simultaneously  results  in  the  equations  for  the  direction 
numbers  a^.  bM,  cjvj,  for  the  line  of  intersection  of  planes  Plane^  and  Planep,  which  are  also 
the  direction  numbers  for  the  principal  axis  through  the  model: 


aM  =  KHLP  ‘  KPLH 

(164) 

bM  =  LHJP  '  LPJH 

(165) 

CM  =  JHKP  ‘  JPKH 

(166) 

The  values  for  JH,  Kh,  are  from  Equations  152  through  155,  and  Jp,  Kp,  Lp  are  from 
Equations  157  through  160,  which  are  defined  in  terms  of  known  values. 
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From  the  direction  numbers  aM,  bM,  cM,  in  Equations  164  through  166,  the  direction 
cosines  of  the  model  principal  axis  are  given  by  the  equations: 


COS  ocjy{  = 


+  /- 


cos  pM  = 


COS  yM 


+/- 


(167) 


(168) 


(169) 


(4)  Find  the  Centroid,  CM,  of  the  Side  View  Area  of  the  Model. 

The  centroid  is  found  as  the  intersection  of  the  two  lines,  lying  in  planes  PlaneH 
and  Planep,  which  are  drawn  from  the  flash  FH  to  the  Hall  shadowgraph  CH  and  from  Fp  to  the 
pit  shadowgraph  Cp.  Point  CM  is  located  by  (XCM,  YCM,  ZcM). 

To  find  the  intersection  of  two  lines  in  3-space,  their  equations  must  be  known  and  solved 
simultaneously.  The  equations  for  each  line  are  written  as  the  equations  for  two  intersecting 
planes,  one  of  which  is  Planepj  for  the  Hall  and  Planep  for  the  Pit.  Equations  for  the  second 
describing  plane  for  each  line  must  be  derived.  The  simplest  plane  in  each  case  will  be  the 
projection  plane  that  is  perpendicular  to  the  range  coordinate  plane  YZ  and  contains  the  line 
F^Cpi  for  the  Hall  and  FpCp  for  the  Pit.  The  general  equation  for  a  line  defined  by  two  points 
Pi,  and  P2,  with  coordinates  Xj,  Yj,  Zj,  and  X2,  Y2,  Z2,  is  : 


and  is  also  defined  as: 


X  -  Xj  Y  -  Yl  Z  -  Z1 
X2  -  xi  Y2  -  Yj  Z2  -  Zj 


x  -  Xj  Y  -  Yj  Z  -  Zj 

~  a  “  b  _  c 


(170) 


(171) 


where  a,b,c  are  the  direction  numbers  for  the  line.  So  we  can  equate: 

a  =X2-Xi 

b  =y2-y, 

c  =  Z2  -  Zj 
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(172) 

(173) 

(174) 


The  equation  for  a  projection  plane,  perpendicular  to  the  YZ  plane  and  passing  through  a 
line,  has  its  X  term  eliminated,  and  the  general  form  reduces  to: 


Y  -  Y, 


(175) 


where  Y  j,  Z[  is  a  point  on  the  line,  and  b,  c  are  direction  numbers  for  the  line  in  the  Y,  Z 
directions.  Rewriting  the  general  form  into  the  more  familiar  equation  for  a  projection  plane 
through  a  line,  with  direction  numbers  a,b,c,  onto  plane  YZ,  we  obtain: 


cY-bZ  +  (bZj  -cYj^O 


(176) 


The  direction  numbers  for  the  lines  FhCh  and  FpCp  are  proportional  to  the  points  on  the 
line  as  in  Equation  171,  so  the  direction  numbers  for  the  Hall  and  Pit  centroid  lines  are: 


aH  =  XCH  -  XFH 
aP  =  XCP  *  XFP 

bH  =  YCH  -  YFH 
bP  =YCP’YFP 

CH  =  ZCH  ZFH 
CP  =  ZCP  •  ZFP 


(177) 

(178) 

(179) 

(180) 

(181) 

(182) 


Substituting  Hall  and  Pit  notations  into  Equation  176  for  b,  c,  we  get  the  two  equations  for  the 
projecting  planes  for  the  Hall  and  Pit  centroid  lines: 


Then  letting: 


cHy  -  hHz  +  (hHZcH  ‘  cHYCH)  =  0 
CpY  -  bpZ  +  (bpZcp  *  CpYcp)  =  0 


Dh  =  bHZcH  *  cHYCH 
Dp  =  bpZcp  -  CpYcp 


we  rewrite  Equations  183  and  184  into  a  simpler  form: 

CpjY  -  bj|Z  +  Dh  =  0 
CpY  -  bpZ  +  Dp  =  0 


(183) 

(184) 


(185) 

(186) 


(187) 

(188) 


where  Equation  187  is  for  the  projection  plane  on  YZ  through  line  FhCh  and  Equation  188  is  for 
the  projection  plane  on  YZ  through  line  FpCp. 

The  intersection  of  lines  FhCh  and  FpCp  can  be  found  by  solving  the  equations  for  two 
planes  forming  fhch  and  one  equation  from  line  FpCp.  This  actually  is  the  intersection  of  line 
FHCH  with  one  plane  describing  line  FpCp.  Therefore,  to  show  that  the  two  lines  do  intersect, 
we  must  also  Find  the  intersection  of  line  fHcH  with  the  other  plane  describing  line  FpCp.  If  the 
coordinates  resulting  from  the  intersection  of  line  FhCh  with  the  two  planes  of  line  FpCp  are  not 
the  same,  we  will  use  for  now  the  average  value  of  the  two  results. 


We  can  solve  for  the  intersection  of  line  fHcH  with  Planep  by  solving  the  equations  for 
line  FhCh  and  Planep  simultaneously. 

Line  is  given  by  Equations  151  and  187: 

JpjX  +  Kp{Y  +  LhZ  +  Mpj  =0  (151) 

chY  -  bHZ  +  Dh  =0  (187) 

Line  FpCp  is  given  by  Equations  156  and  188: 

JPX  +  KPY  +  LPZ  +  MP  =0  (156) 

CpY-bpZ  +  Dp  =0  (188) 

The  general  solution  to  three  equations  with  three  unknowns,  X,Y,Z  from  the  form: 


is  given  by: 


A1X  +  B1Y  +  C1Z  =  Dj 
A2X  +  B2Y  +  C2Z  =  D2 
A3X  +  B3Y  +  C3Z  —  D3 


Y  AX  v  AY  „  _  AZ 
A"A’Y“A,Z'“A 


A  ’  A’  A 
where  the  determinants  for  the  solutions  of  A,  AX,  AY,  AZ  are: 


(189) 

(190) 

(191) 

(192) 


At 

Bi 

Ci 

Di 

Bi 

c, 

a2 

b2 

c2 

(193) 

AX  =  D2 

b2 

c2 

a3 

«3 

c3 

d3 

b3 

c3 

(194) 


Ai 

Di 

Cl 

Ai 

Bi 

D, 

a2 

d2 

C2 

(195) 

AZ  =  A2 

b2 

d2 

a3 

D3 

C3 

a3 

b3 

d3 

(196) 


LA.  K  * 


Solving  Equations  151,  187,  and  156,  we  get  the  A  values,  denoted  by  subscript  155,  for 
the  intersection  of  line  F^jC^  with  Planep  using  the  general  solution  to  three  equations  with  three 
unknowns  described  above: 


A156  = 

Jh 

0 

JP 

kH  lH 

CH  'bH 
KP  LP 

(197) 

^156  = 

-Mh  kh  lh 

-E>h  cH  -bH 
-Mp  Kp  Lp 

(198) 

AY  !56  = 

Jh 

0 

Jp 

-Mh  lh 
-Dh  "bH 
-Mp  Lp 

(199) 

AZ156  = 

jh  kh  -mh 
0  ch  -Dh 
Jp  Kp  -Mp 

(200) 

Likewise,  solving  Equations  151,  187,  and  188,  we  get  the  A  values,  denoted  by  subscript 
jgg,  for  the  line  FHCH  intersection  with  the  projection  plane  onto  YZ  through  line  FpCp: 

A 188  = 

Jh 

0 

0 

Kh  lh 

CH  *bH 

Cp  -bp 

(201) 

AX188  = 

-mh  Kh  Lh 

-Dh  ch  _bH 
-Dp  Cp  -bp 

(202) 

AY 188  = 

jh 

0 

0 

-mh  lh 

-Dh  -bH 
-Dp  bp 

(203) 

AZl8L  = 

JH  kh  -mh 
0  ch  -Dh 

0  Cp  -Dp 

(204) 

Using  the  solutions  for  the  A  values  from  Equations  197  through  200  for  the  intersection  of  line 
FhCh  and  Planep  (A  156  subscripts)  and  the  solutions  for  the  A  values  from  Equations  201 
through  204  for  the  intersection  of  line  FHCH  and  projection  plane  onto  YZ  through  line  FpCp 
(A  jgg  subscripts),  we  can  write  the  equations  for  the  coordinates  of  the  centroid  of  the  model: 

on  Planep 

Cm  Plane.  YZ-Projegiion 

XCM156  = 

AX>56 

A156 

(205) 

XCM188 

AX188 

A188 

(208) 

YCM156  = 

AY.56 

A156 

(206) 

YCM188 

AY188 

A188 

(209) 

ZCM156 

AZ156 

A156 

(207) 

ZCM188 

AZ188 

A188 

(210) 

If  the  values  from  Equations  205  through  207  are  different  from  those  of  Equations  208 
through  210,  then  the  values  to  use  for  the  location  of  the  centroid  of  the  model  should  be  the 
average  of  the  values.  (In  reality  the  centroid  point  should  be  half  the  length  of  a  perpendicular 
drawn  between  the  closest  points  of  the  two  lines,  but  we  will  use  the  average  for  now.) 

So  the  average  values  for  the  coordinates  of  the  model  centroid  Cps/f  are: 


*CM  156  +  ^CM188 

XCM  = - 2 -  (211^ 

Y  +  Y 
CM  156  CM188 

YCM  =  - 3 -  (212) 

2 CM  156  +  ^TM188 

Zcm  = - 2 -  (213) 

The  coordinate  locations  for  the  center  of  mass,  CG,  of  the  model  are: 

XCG  =  XCM  +  Ax  (214) 

YCG  =  YcM  +  Ay  (215) 

ZcG  =  Zcm  +  Az  (216) 

The  Ax,  Ay,  Az  values  are  related  to  the  distance,  dCG,  between  the  centroid  and  the  center 
of  mass  along  the  model  axis,  which  is: 

dec  =  (Mg  -  Mco)  (2 1 7) 


where  Me  and  Mcg  are  known  values  obtained  from  the  model.  From  simple  trigonometry,  the 
equations  for  Ax,  Ay,  Az  are: 


Ax  =  dGG  COS  aM 

(218) 

Ay  =  deG  COS 

(219) 

Az  =  deG  COS  Ym 

(220) 

Substituting  Equations  217  through  220  into  Equations  214  through  216, 
coordinate  locations  for  the  center  of  mass,  CG,  of  the  model: 

we  get  the 

Xcg  =  Xcm  +  dec  COS  aM 

(221) 

Ycg  =  yCM  +  dcG  COS  Pm 

(222) 

Zcg  =  Zcm  +  dcG  COS  Ym 

(223) 
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7.  COMMENTS 


The  solutions  for  the  centroid,  center  of  mass,  and  direction  cosines  of  the  principal  axis 
through  the  model  were  derived  using  known  constants  for  the  locations  of  the  reference  beads  in 
both  the  Hall  and  Pit  views.  The  location  of  Bead  Bjh  does  not  necessarily  have  to  be  the 
coordinate  location  (0,0,0),  so  the  equations  are  general  enough  to  be  used  for  any  station  along 
the  range. 

The  location  for  the  CG  of  the  model  is  based  on  the  assumption  that  the  CG  and  the 
centroid  of  the  projection  of  the  model  lie  on  the  same  major  axis.  If  they  do  not,  Equations  22 1 
through  223  must  be  modified. 


The  spatial  resolution  factors  determined  for  the  Hall  and  Pit  views  are  only  the  horizontal 
resolutions,  determined  by  the  measurement  and  known  distance  between  beads  1  and  2  in  each 
view.  If  the  vertical  resolution  is  to  be  used  also,  the  camera  measurement  of  the  distance 
between  bead  1  and  bead  3  must  be  used.  However,  Eglin  has  found  that  the  vertical  and 
horizontal  resolutions  are  almost  always  equal. 
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SECTION  VII 


PROOF-OF-PRINCIPLE  SYSTEM 


The  proof-of-principle  system  for  high-speed  electronic  imaging  was  delivered  to  the 
Aeroballistics  Research  Facility  of  the  Eglin  Air  Force  Base  on  May  18,  1987. 

The  purpose  of  the  system  was  to  demonstrate  the  feasibility  of  either  replacing  or 
augmenting  the  current  film  camera  system  presently  operating  in  the  Research  Facility  with  solid 
state  video  cameras.  This  demonstration  was  to  show  the  measurement  resolution  capability  for 
locating  projectiles  as  they  fly  down  the  range  and  also  to  demonstrate  proper  sensitivity  and 
operation  in  situ  at  the  range  using  existing  range  timing,  controls,  and  spark-gap  illumination. 

Another  purpose  of  the  demonstration  was  to  show  the  effects  of  special  purpose 
algorithms  and  image  processing  on  improving  the  accuracy  of  the  measured  projectile  location. 
The  goal  was  to  show  the  potential  of  matching  the  accuracy  of  the  film  system. 

Primary  reasons  to  replace  the  film  system  in  the  range  with  electronic  imaging  is  to  reduce 
the  time  to  get  test  data  from  6  to  8  weeks  to  a  few  hours,  eliminate  the  installation's  dependency 
on  single-source  film,  and  to  improve  repeatability  of  measurements. 

I.  CONFIGURATION 

The  study  detailed  in  Section  V  has  shown  that  1024v  x  512h  resolution  electronic  imaging 
cameras  have  sufficient  accuracy,  when  used  with  centroid  and  second-moment-of-inertia 
algorithms,  to  yield  measurement  results  comparable  to  the  film  system.  Based  upon  that,  the 
proof-of-principle  system  simulated  a  1024v  x  512h  camera  system  by  using  zoom  lenses  on  512 
x  5 12  resolution  cameras  to  reduce  the  field  of  view  to  achieve  equal  spatial  resolution. 

The  1024v  x  512j,  cameras  had  to  be  simulated,  because  a  camera  with  the  right  resolution 
and  price  was  not  currently  available.  The  proof-of-principle  system  was  then  configured  as 
follows: 

a .  Two  GE  Model  2250  charge-injection  device  (CED)  Cameras  with  12.5  to  75-mm 
zoom  lenses,  5 12v  x  5 12^  resolution,  single-frame  integrate-mode  operation. 

b.  Two  Poynting  Products  Model  M5250-PC  controllers  for  the  GE  cameras, 
interfaced  to  an  IBM  PC/AT  computer. 

c .  IBM  PC/AT  computer  with  1  -mb  memory,  running  PC-DOS  vs  3. 1 ,  with  a 
monochrome  monitor. 

d  Imaging  Technology  Model  FG- 100- AT- 1024  image  processing  card,  with  1  mb 
memory.  The  PC  card  was  installed  in  the  IBM  AT. 

e.  Sony  19"  color  monitor  with  51 2h  x  490v  resolution.  The  color  monitor  was 
connected  to  the  FG-100  PC  card. 
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f .  Wemer-Frei  "Imagelab"  image  processing  software  (Reference  7)  for  the  FG-100. 
This  includes  a  subroutine  library  "Toolbox"  with  subroutines  usable  from  Fortran 
programs. 

g  .  Microsoft  Mouse  and  Mouse  driver  software.  This  is  a  PC  card  installed  in  the 
IBM  AT,  along  with  the  driver  software.  The  Mouse  is  used  by  the  Wemer-Frei 
"Imagelab"  operating  system  to  select  menu  commands  or  areas  in  an  image  for 
processing. 

h .  Infrared  Filters  installed  in  the  zoom  lenses  to  filter  out  the  IR  screen  light  used  at 
the  facility.  The  filter  cutoff  is  about  700  nm. 

2.  ASSEMBLY  AT  LOS  ALAMOS 

The  system  was  assembled  and  tested  at  Los  Alamos  prior  to  delivery  to  Eglin.  Additional 
software  was  also  written  to  provide  some  image  processing  capability  not  provided  by  Wemer- 
Frei  .  The  two  algorithms  (centroid  and  second-moment-of-inertia)  were  also  coded,  along  with 
the  equations  for  locating  an  object  in  3-space  from  two  orthogonal  point  projections.  Software 
to  operate  the  cameras  was  also  written,  and  menus  for  operation  were  embedded  in  the 
"Imagelab"  operating  system. 

The  system  and  software  were  tested  by  operating  the  cameras,  processing  images,  and  by 
calculating  the  centroids  and  second-moments-of-inertia  of  several  objects  of  known  size  and 
orientation  placed  directly  in  memory  (not  imaged  by  a  camera).  The  algorithms  and  software  all 
seemed  to  work  with  very  high  accuracy. 

3  .  INSTALLATION  AT  EGLIN 

On  May  18,  1987,  the  system  was  delivered  to  Eglin.  It  was  unpacked,  assembled,  and 
preliminary  tests  found  it  to  be  in  good  operating  order.  Mr.  R.  J.  Kelley,  manager  of  the 
research  facility,  allocated  space  and  locations  for  the  equipment,  and  assigned  John  Krieger, 
technician,  to  us  to  help  in  the  installation.  This  included  running  cables,  making  connections, 
fabrication  of  camera  and  mirror  mounts,  and  connecting  the  camera  system  to  the  range  firing 
control  system.  Figure  19  shows  the  IBM  AT  computer  and  the  Sony  color  monitor  located  in 
the  control  room  at  the  Aeroballistics  Research  Facility.  Cables  were  run  overhead  in  cable  trays 
to  the  cameras  in  the  firing  range  and  to  the  control  console  in  the  same  control  room. 


The  Poynting  Products  camera  controllers  were  located  in  the  Public  Address  amplifier 
rack,  behind  the  computer.  It  should  be  mentioned  that  the  cables  between  the  cameras  and 
controllers  were  100  feet  long,  a  factor  which  resulted  in  some  degradation  of  the  images.  The 
full-range  implementation  will  have  the  image  buffers  in  the  camera  itself,  eliminating  this  long 
cable  run  and  its  resultant  noise  and  loss  of  resolution. 

Figure  20  is  a  photograph  of  the  Hall  camera  installation  at  Station  1.  The  film  camera  is  at 
the  top,  the  spark  gap  is  the  circle  in  the  center,  and  the  IR  screen  is  at  the  bottom  right  of  the 
picture.  The  GE  camera  is  at  the  lower  center  and  is  mounted  closer  to  the  spark  gap  than  the 
film  camera  to  compensate  for  the  lower  CID  sensitivity.  Tests  have  shown  that  the  CCD  camera 
is  about  three  lens  f-stops  more  sensitive  than  the  GE  CED  cameras.  This  picture  shows  that  the 
electronic  imaging  camera  can  be  installed  alongside  the  existing  range  cameras  and  equipment 
without  disturbing  their  installation.  This  allows  the  film  system  to  remain  intact  and  operational 
while  the  electronic  imaging  system  is  installed. 

Figure  21  shows  the  Pit  camera  installation,  which  was  slightly  different  from  the  Hall 
installation,  in  that  a  turning  mirror  was  used  because  of  the  different  mount  configuration  for  the 
Pit  camera.  This  causes  the  image  to  be  backwards,  but  that  was  reversed  by  the  computer. 

Part  of  the  spark  gap  can  be  seen  between  the  film  camera  and  the  mirror.  The  Pit  does  not 
have  an  IR  screen,  because  its  spark  gap  is  triggered  from  the  corresponding  Hall  IR  screen. 
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4  .  CONTROLS  AND  OPERATION 

Controls  for  the  electronic  imaging  system  were  simple.  The  operator  selected  the  ARM 
CAMERAS  command  from  the  menu  on  the  monitor,  and  the  computer  set  the  cameras  into  the 
integrate  mode  of  acquiring  a  single  frame.  This  was  done  twenty  seconds  prior  to  firing  time. 
Two  control  signals  from  the  range  control  panel  were  used  to  start  and  stop  the  camera 
integration  period,  similar  to  opening  and  closing  a  film  camera  shutter.  These  were  the  control 
panel  Lockout  On  and  Lockout  Off  signals,  with  their  actuation  timer  setting  from  thumbwheel 
switches  on  the  panel.  To  allow  for  electro-mechanical  delays  in  the  gun  firing  mechanism,  the 
electrical  signal  firing  pulse  occurs  1.000  second  before  T=0.  The  Lockout  On  signal  was  also 
set  at  -1.000  second,  and  the  Lockout  Off  signal  at  -0.700  second.  As  soon  as  the  shot  was 
fired,  the  computer  was  then  accessed,  and  an  image  put  on  the  monitor  within  30  seconds. 

5 .  LIVE  TESTS 

During  the  extent  of  the  installation  and  testing  period,  a  number  of  images  of  live  models 
fired  down  the  range  were  captured  and  stored.  These  were  often  used  by  the  range  personnel 
for  "quick  look"  information  to  determine  if  the  shot  was  good  or  not.  Unfortunately,  we  did  not 
have  enough  time  to  test  our  algorithms  and  software  for  locating  the  model  in  3-space,  because 
all  our  time  was  taken  in  other  testing  and  image  processing.  Completing  these  tests  is  a  high 
priority  and  should  be  conducted  as  soon  as  sufficient  funding  is  available,  to  prove  the 
algorithms  and  accuracy  of  the  location  calculations.  The  images  were  stored  on  discs  for  future 
reference. 

6.  1024v  x  512h  CCD  CAMERA  SYSTEM 

In  conjunction  with  the  studies  and  algorithm  development  that  Los  Alamos  has  done  for 
Eglin,  we  invited  a  vendor  to  demonstrate  a  prototype  1024v  x  5 12]-,  CCD  camera  at  Eglin  during 

the  testing  of  the  proof-of-principle.  We  had  initially  made  contact  with  the  vendor  through 
Megavision.  Because  the  Megavision  Image  Processor  was  already  interfaced  to  the  CCD 
camera,  we  also  invited  Megavision  to  capture  some  images  on  the  Eglin  test  range  so  that  we 
might  compare  them  with  the  proof-of-principle  images  and  also  to  verify  that  the  CCD  cameras 
were  sensitive  enough  to  work  with  the  Eglin  spark  gaps. 

It  took  several  days  to  get  the  Megavision/CCD  Camera  system  operating  with  the  range 
control  signals,  but  we  did  manage  to  capture  one  test  image  of  a  projectile  fired  in  the  range. 

We  took  many  other  images  of  static  projectiles  for  calibration  purposes  and  later  transferred 
those  images  to  the  IBM  AT  system,  where  we  performed  some  calculations  on  one  image.  The 
Megavision  system  provided  an  easily  adaptable  interface  to  the  CCD  camera  and  also 
demonstrated  image  processing  capability  that  would  provide  the  necessary  functions  for  the 
full-range  instrumentation.  It  would  be  to  Eglin 's  advantage  to  do  all  further  algorithm  and 
image  processing  development  on  the  image  processor  that  will  be  used  on  the  full  range  system 
so  that  all  software  developed  will  be  directly  usable  and  not  wasted.  A  small  system  like  the 
IBM  AT  does  not  have  sufficient  capabilities  for  the  range.  The  Megavision  system  is  interfaced 
to  the  CCD  camera,  has  the  necessary  capabilities,  and  would  be  a  good  choice. 

Even  though  the  CCD  camera  was  a  prototype  and  did  not  contain  a  controller  or  frame 
buffer,  it  did  show  very  good  performance  that  appears  to  meet  the  range-imaging  requirements. 
Characterization  and  more  testing  of  the  CCD  camera  should  be  done  next  year  to  support  this 
early  work  before  committing  to  it  for  the  full  range  instrumentation. 
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7 .  ALGORITHM  TESTS 

Prior  to  delivery  of  the  proof-of-principle  system  to  Eglin,  the  centroid  and  second- 
moment-of-inertia  algorithms  were  coded  and  tested  at  Los  Alamos. 

a.  Centroid 

Because  the  bead  locations  are  very  important  in  any  measurement  and  because  they 
are  also  very  small  in  size,  we  simulated  a  bead  in  image  memory  and  used  that  to  test  the 
centroid  algorithm.  The  detailed  Eglin  bead  is  shown  in  Figure  6. 

For  a  1024v  x  5 12^  pixel  camera,  the  spatial  resolution  for  a  96v”  x  48h"  field  of  view 
(the  reflecting  screen)  is  the  same  for  both  the  horizontal  and  vertical  directions,  so  we  can 
calculate  it  using  the  horizontal  dimensions  as: 


Rs  " 


Field  of  view  dimension 
Number  of  pixels 


48" 


0.09375" 


512  pixel  pixel 


(224) 


Using  Rs  for  both  pixel  dimensions,  a  pixel  area  is  then  R2  =  8.789  x  10~3  _£uare  mc^s 

s  square  pixel 

We  will  also  assume  a  viewing  situation  with  the  camera  that  should  be  the  worst  case,  i.e.,  most 
pixels  in  the  camera  are  only  partially  covered  by  the  image,  as  shown  in  Figure  22.  We  will, 
however,  arrange  the  bead  such  that  we  can  calculate  the  sub-pixel  areas  conveniently. 


Figure  22.  Eglin  Bead  Image  on  CCD 

The  pixels  covered  by  the  bead  are  numbered  1  through  12.  The  calculated  sub-pixel  areas 
of  the  bead  image  are  listed  in  Table  3,  along  with  the  image  processor  digitized  values  rounded 
to  8-bit  integers  (the  8-bit  digitizer  values  are  0  for  zero  area  and  254  for  an  area  of  8.789  x  10'-' 
square  inches,  or  one  square  pixel). 
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TABLE  3.  BEAD  PIXEL  AREAS 


PIXEL 

AREA 

VALLE 

1 

.002180 

63 

2 

.003176 

92 

3 

.001512 

44 

4 

.000443 

13 

5 

.003000 

87 

6 

.008789 

254 

7 

.008789 

254 

8 

.002578 

75 

9 

.002180 

63 

10 

.003175 

92 

11 

.001512 

44 

12 

.000443 

13 

The  digitizer  values  calculated  for  the  bead  were  then  entered  into  the  image  processor 
memory,  resulting  in  an  image  shown  in  Figure  23  (magnified),  where  pixel  values  of  254  are 
white  and  values  of  0  are  black. 
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Figure  23.  Magnified  Bead  Image  on  Monitor 


The  centroid  of  the  bead  has  been  calculated  manually  as  being  along  the  central  horizontal 
axis  and  0.1074"  from  the  left  edge  (top  hat  section).  This  was  done  by  calculating  the  centroids 
of  the  two  distinct  rectangular  areas  of  the  bead,  summing  the  individual  moments  of  the  two 
rectangular  areas,  and  then  dividing  by  the  total  area,  or 


X 


c 


Sum  of  Area  Moments  0.004056559 
Total  Area  “  0.037778 


0.1073789 


(225) 


Rounding  Xc  to  four  places  results  in 

Xc  =  0.1074 


(226) 


In  the  image  processor,  the  center  of  pixel  1  is  location  0,0  for  x  and  y,  where  y  is  positive 
down  and  x  is  positive  to  the  right.  The  location  of  the  left  edge  of  the  bead  from  the  center  of 
pixel  5  is  then  calculated  as 


v  _  0.09375" 

A(leftedge)  ~  2 


0.032" 


0.014875" 


(227) 
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Then,  because  the  centroid  of  the  bead  is  0.1073789"  from  the  left  edge,  the  centroid  of  the  bead 
relative  to  the  center  of  pixel  5  (location  0,1)  is 

Xc  (inches)  =  0.014875  +  0.1073789  =  0.1222539"  (228) 

Changing  this  dimension  from  inches  to  pixels  (which  is  what  the  IBM  image  processor  uses): 


0 1222539" 

(theoretical)  Xc  (pixels)  =  o.0937?7pixeI  =1-30404  Plxels 

This  is  the  theoretical  location  of  the  X  dimension  at  the  bead  centroid.  The  Yc  (pixels)  ’s 
calculated  similarly  as: 


(229) 


(theoretical)  Yc  (pixels)  =  1-0  pixels 


(230) 


The  image  processor  was  then  allowed  to  calculate  the  centroid  of  the  simulated  bead  image 
so  we  could  determine  the  highest  accuracy  of  the  algorithm.  This  is  only  a  theoretical 
comparison,  because  the  bead  image  used  contains  absolutely  no  noise  (the  presence  of  which 
will  reduce  the  accuracy).  In  the  real  world  we  do  have  noise,  so  measurements  were  also  done 
at  Eglin  using  real  images.  But  it  is  useful  to  know  the  ultimate  theoretical  precision  of  the 
programmed  algorithm  in  the  system. 

The  values  calculated  by  the  IBM  AT  image  processor  using  the  programmed  centroid 
algorithm  are  tabulated  in  Table  4,  with  the  theoretical  values. 

TABLE  4.  CALCULATED  BEAD  AREA  CENTROIDS 


Theoretical  Xc 
Calculated  Xc 
Xc  Error 
Theoretical  Yc 
Calculated  Yc 
Yc  Error 


1.30426  pixel 
1.3015  pixel 
0.00276  pixel 
1 .0  pixel 
1 .0  pixel 
0  pixel 


It  may  be  seen  that  the  largest  error  was  in  the  X  dimension  =  0.00276  pixel.  Because  the 
spatial  resolution  in  the  image  is  0.09375  "/pixel,  the  error  is  0.00276  x  0.09375,  or  Xc(Error) 

=  0.000259",  which  is  extremely  good  relative  to  the  Aeroballistics  range  requirements.  The 
claimed  best  case  error  for  the  presently  used  film  camera  system  on  the  range  is  0.03  cm 
(0.012")  for  the  projectile  position  (Reference  8).  This  translates  to  about  0.024"  best  case 
shadowgraph  error  for  the  film  system  by  assuming  the  projectile  in  the  center  of  the  range  and 
using  simple  geometry.  This  agrees  closely  with  a  calculated  film  error  of  about  0.023" 
determined  by  the  Los  Alamos  study  done  in  1985  (Reference  5). 

b.  Angle 

A  process  similar  to  the  centroid  comparison  was  done  to  test  the  second-moment-of- 
inertia  (angle)  algorithm  in  the  computer.  A  small  rectangular  shape  approximately  1 .875"  long  \ 
0.46875"  wide  was  chosen.  The  rectangular  shape  was  oriented  at  an  angle  corresponding  to  a 
rise  of  Y  =  1,  X  =  5  or 

^(theoretical)  =  TAN  ^  J  =  11.30993247  (_31) 
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The  sub-pixel  values  could  then  be  calculated  precisely  because  the  rectangle  image  was 
located  advantageously  on  the  imager.  The  IBM  AT  image  processor  second-moment-of-inertia 
algorithm  software  calculated  an  angle  as 

9  (calculated)  =  11-30992°  (232) 

The  difference  between  the  theoretical  and  calculated  angles  is  in  error  of  only  0.00001  °.  This 
again  is  for  the  case  with  no  noise.  But  because  the  angle  algorithm  also  calculates  and  uses  the 
centroid,  we  see  that  the  small  angle  error  is  very  good  indeed. 

8.  MEASUREMENTS  AT  EGLIN 

The  only  measurements  taken  at  Eglin,  limited  by  the  short  amount  time  we  had  to  do  other 
testing,  imaging,  and  installation,  were  centroids  of  beads  and  spatial  locations  of  reference 
objects.  These  reference  objects  were  thumbtacks  mounted  on  the  screen  itself.  The  tack 
locations  relative  to  the  beads  are  shown  in  Figure  24. 


Figure  24.  Bead  and  Tack  Locations 

The  diameters  of  the  two  tacks,  designated  as  points  5  (tack  2)  and  6  (tack  1)  in  the  image,  were 
measured  as: 

DS  =  0.3697”  and  D6  =  0.3760"  (233) 

Figure  25  is  a  photograph  of  the  Sony  monitor  showing  an  image  of  the  screen  taken  with 
the  CCD  camera.  The  two  bullet  shapes  are  a  projectile  and  its  shadowgraph.  Tacks  1  (point  6) 
and  2  (point  5)  and  beads  1  and  2  are  indicated.  This  image  is  one  of  two  used  for  this 
measurement. 

The  second  image  was  taken  with  the  GE  camera  and  is  similar  to  the  first,  but  without  the 
projectile.  We  did  not  photograph  that  image. 
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a.  Plan 

Ideally,  we  wanted  to  very  accurately  measure  the  distance  between  the  two 
thumbtacks  to  use  as  a  standard  for  comparison  with  the  same  measurement  obtained  by  the 
image  processing  algorithms  run  on  the  EBM  AT  and  by  the  film  cameras.  The  plan  was  to: 

( 1 )  Measure  the  distance  between  the  tacks  using  a  tape  measure.  This  is  a  coarse 

measurement,  but  it  yields  ~  (0.0625")  accuracy. 

(2 )  Photograph  the  same  screen  with  the  existing  film  camera  and  measure  the 
distance  between  the  tacks  using  the  normal  film  reading  procedure  of  the  Math  Laboratory  at 
Eglin.  We  are  not  sure  of  the  accuracy  of  this  measurement,  but  it  should  have  a  probable  error 
of  about  +/-  0.023",  under  ideal  conditions,  because  that  number  was  determined  by  an 
independent  study  (Reference  5)  by  Los  Alamos  in  1985  by  measuring  the  modulation  transfer 
function  (MTF)  of  the  film  camera  system.  Of  course,  different  cameras,  lighting,  screen 
conditions,  and  location  of  the  bead  in  the  image  all  affect  this  accuracy.  But  for  now,  we  will 
use  that  measurement  of  0.023". 

(3 )  Measure  the  distance  between  the  tacks  with  the  calibration  transit  used  by  the 
Aeroballistics  Research  Facility  in  making  range  measurements  of  actual  bead  locations.  This 
system  is  probably  accurate  to  about  +/-  0.005". 

(4 )  Measure  the  distance  between  the  tacks  with  the  GE  and  CCD  cameras  and 
using  the  centroid  algorithm. 

b .  Film  Measurement  of  Tacks 

A  photograph  of  the  screen  was  taken  and  sent  to  the  Math  Laboratory  at  Eglin  for 
reading  in  the  normal  fashion,  except  points  5  and  6  were  added,  as  shown  in  Figure  25. 

Locations  of  points  5  (tack  2)  and  6  (tack  1)  are  the  inside  edges,  just  as  the  beads  are 
measured.  Results  of  their  measurements,  in  microns,  are  tabulated  in  Table  5. 


TABLE  5.  MATH  LABORATORY  FILM  MEASUREMENTS 


Of  course,  for  statistically  more  reliable  data  ,  this  measurement  should  be  repeated  10  to 
15  times.  But,  in  our  case,  it  was  read  only  once  during  the  time  we  did  this  testing. 

We  know  that  the  inside  edges  of  the  beads  are  set  36.000”  apart  when  they  are  assembled 
on  the  Kevlar  line.  Any  error  in  that  assembly  will  cause  an  absolute  error  in  measurements  in 
the  range.  But  the  beads  are  assembled  very  carefully  on  pre-stretched  line  using  a  machined 
fixture,  so,  for  now,  we  will  assume  no  error  in  this  measurement.  Measuring  the  distance 
between  inside  edges  of  the  beads  on  the  film,  and  using  36.000"  as  the  actual  distance,  allows 
us  to  establish  the  spatial  resolution  factor  for  other  film  measurements. 

Because  the  beads  and  tacks  in  the  image  may  be  tilted  with  respect  to  the  camera,  we  can 
use  simple  geometry  to  find  the  straight-line  distance  between  them.  For  the  beads,  is  the 
straight  line  distance  shown  in  Figure  26. 

Bead  1 


Bead  2 


Figure  26.  Bead  Location  Geometry 


where 


ab  =  YB2-  YbI 

bp  =  XB2-Xsi 


rb- A/ab  +  bb 


So,  from  the  film  measurements,  we  find 


ab  =  23197 -  22699  =  -498 
bb  =  72768  -  25395  =  47373 


rb  =  -v  /(- 498)  +  (47373)“  =  47375.62 


The  spatial  resolution  for  the  film,  then,  is: 


„  rb  47375.62  micron 

Rfllm  "  36.000"  _  36.000  -1315  989444  inch 


t  *  »>  ,V’»  ,V, 


V.'/  v.V-V.v.Y./.v.y  j-.j-  .-  w  w 


Also  for  the  tacks,  using  the  same  geometrical  calculation  for  rj: 

aT  =  63225  -  63947  =  -722 
br  =  71997 -  25318  =  46,679 


and 


rT 


122)2  +  (466?9)2 


46684.58 


(241) 

(242) 

(243) 

(244) 


Then  applying  Rfi[m  from  Equation  240  to  the  measurement  rT  to  find  the  distance  between 
the  tacks  measured  on  film,  we  get 

^  rT 

DT(film)  =  -  = 

Rfilm 

DT(film)  =  35.475  " 

c.  Electronic  Imaging  Measurement 

Measuring  the  tack  and  bead  positions  using  the  centroid  algorithm  is  not  as 
straightforward  as  it  may  seem,  because  of  the  very  important  influence  of  noise  in  the  image.  In 
this  case,  noise  is  anything  other  than  the  object  to  be  measured.  Thus,  the  Kevlar  line, 
projectile,  handprints,  and  random  noise  (both  optical  and  electrical)  are  all  "noise."  The 
procedure,  then,  is  to  eliminate  all  extraneous  information  in  the  image,  taking  care  not  to  destroy 
any  measurement  information  in  the  object.  Any  image  processing  filters  or  thresholding  applied 
to  the  image  must  be  such  as  to  comply  with  the  stated  procedure.  Some  of  the  experimenting  at 
Eglin  was  an  effort  to  derive  some  filtering  methods  to  eliminate  all  the  background  (set  to  zero 
value)  and  to  leave  the  object  (tack  or  bead)  intact.  Some  procedures  were  derived  by  trial  and 
error  using  knowledge  of  what  affects  measurement  information,  but  much  image  processing 
development  for  our  measurements  needs  to  be  done.  This  seems  to  be  a  relatively  new,  or  at 
least  unpublished,  area  of  image  processing;  that  is,  most  image  processing  has  to  do  with 
improving  how  an  image  looks  to  the  human  eye,  but  not  in  improving  measurement  capability 
within  an  image. 

Taking  some  simple  tries  at  processing  an  image  for  eliminating  the  background,  two 
different  images  acquired  by  two  different  cameras  and  digitizing  systems  were  operated  on  to 
measure  the  distance  between  the  tacks. 


46684.58 

1'3 15.989375' 


of  Tacks 


(245) 

(246) 
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(1)  GE  Camera  Image 

The  centroid  data  measured  on  the  GE  camera  image  is  shown  in  Table  6,  in 

pixel  units  as  : 

TABLE  6.  CENTROID  MEASUREMENTS  WITH  GE  CAMERA 


Calculating  the  image  spatial  resolution  using  the  same  geometry  as  for  the  film  (Figure  26),  we 
get: 


ab(GE)  =  446.3474  -  447.7971  =  -  1.4497 

(247) 

bb(GE)  =  36.4042  -  459.9247  =  -  423.5205 

(248) 

i-b(GE)  =  ^aj;  +  b£  =  ^(-1.4497)Z+  (-423.5205)2 

(249) 

rb(GE)  =  423.52298 

(250) 

However,  r^g)  for  the  electronic  imaging  cameras  is  the  distance  measured  between 
centroids  of  the  beads  (see  Figure  6),  which  should  be  36.000"  +  2(0.1074")  =  36.2148".  So 
the  spatial  resolution  on  the  GE  camera  image  is: 

n  423.52298  ,,  pixel 

R<ge>  -  Mil  "  1-69474859  ^ 

(251) 

The  tack  measurements  are  : 

aT(GE)  =  37.3495  -  86.3240  =  1.0255 

(252) 

bT(GE)  =  47.4939  -  466.4751  =  -  418.9812 

(253) 

rT(GE)  -  +  b  T 

(254) 

rT(GE)  =  4 1 8.982455 1 

(255) 

Applying  the  spatial  resolution  of  the  image  to  the  distance  between  centroids  of  the  tacks,  we  get 

418.9824551  pixel  „  „„ 

dTC(GE)  = - =  35.826;>  (256) 

1 169474859 


Subtracting  half  of  the  two  tack  diameters  will  give  us  the  distance  between  the  inside  edges  of 
tacks,  we  get 


DT(GE)  =  35.3265  -  ^ 
dT(GE)  =  35.454" 


0.3697  0.3760 


=  35.4537" 


(257) 

(258) 


units,  as  : 


(2)  CCD  Camera  Image 

The  centroid  data  measured  on  the  CCD  camera  image  is  shown  in  Table  7,  in  pixel 

TABLE  7.  CENTROID  MEASUREMENTS  WITH  CCD  CAMERA 


Object 

X  (pixels) 

Y  (Pixels) 

Bead  1 

459.7842 

447.1870 

Bead  2 

36.8718 

446.0997 

Tack  1 

446.8079 

86.271) 

Tack  2 

47.2950 

87.7229 

Calculating  the  image  spatial  resolution  using  the  same  geometry  as  for  the  film  (Figure  26),  we 
find  : 


ab(CCD)  =  446.0997  -  447.1870  =  -  1.0873 
bb(CCD)  =  36.8718  -  459.7842  =  -  422.9124 

rb(CCD)  =  ^ab  +  bb  =  V(-1.0873)2  +  (422.9124)* 
rb(CCD)  =  422.9137976 

The  spatial  resolution  on  the  CCD  camera  image  is 


(259) 

(260) 

(261) 

(262) 


„  422.913796 

R(CCD)  ~  36.2148 


11.67792719  fe- 
mch 


IS; 
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The  tack  measurements  are: 


aT(CCD)  =  87.7229  -  86.27 1 1  =  1 .45 1 8 
bT(ccD)  =  47.2950  -  465.8079  =  -  418.5129 


(264 


(265) 


rT(CCD)  =  '\jaT  +  bT 

(266) 

rT(CCD)  =  418.515412  pixels 

(267) 

we  get : 


Applying  the  spatial  resolution  of  the’ image  to  the  distance  between  centroids  of  the  tacks. 

(268) 


DtC.CCD,  -  418-51?412S  -  35.8382" 

1 1.67792717  UiS- 
inch 


f 


Subtracting  half  of  the  two  tack  diameters  will  give  us  the  distance  between  the  inside  edges  of 
the  tacks,  we  get 


r-x  „  o-xv-x  0.3697  0.3760  Ae.c^ 

Dt(CCD)  =  35.8382  - 2 - 2 —  =  35.4653  (269) 

DT(CCD)  =  35.465”  (270) 


(3)  Tabulation  of  Measurement  Results 

Shown  in  Table  8  are  results  of  the  measurement  of  the  spacing  between  the 
inside  edges  of  the  tacks  as: 


TABLE  8.  MEASUREMENT  RESULTS 


Method 

Tack  Measurement  1 

Tape 

35.406  " 

Film 

35.476  " 

GECID 

35.454  " 

CCD 

35.465  " 

Transit 

35.380  " 

9.  SUMMARY 


In  summary,  the  electronic  imaging  camera  measurements  resulted  in  differences  of  0.022" 
and  0.0 11"  from  the  film  system  measurement.  The  error  for  the  film  measurement,  using  the 
transit  measurement  as  the  standard,  was  0.096".  The  errors  for  the  GE  and  CCD  cameras  were 
0.074"  and  0.085",  respectively,  with  the  electronic  shadowgraph  measurements  from  both 
cameras  closer  to  the  tape  and  transit  measurements  than  was  the  film  measurement.  Not  too 
much  should  be  made  of  these  results,  because  they  are  based  on  so  few  measurements  and  they 
were  fairly  coarse.  Two  electronic  imaging  measurements,  with  pixel  resolutions  of  about 
0.085"  agree  within  0.011"  and  they  were  only  0.011"  and  0.022"  different  from  the  film 
measurement,  which  contains  some  error.  Since  the  film  measurements  are  so  much  in  error 
from  the  transit  measurements,  something  very  basic  could  be  wrong  with  the  system  or  the 
measuring  devices.  However,  these  results  for  the  electronic  cameras  look  very  encouraging, 
and  certainly  are  within  the  realm  of  accuracy  of  the  present  film  camera  system,  with  a  distinct 
possibility  of  being  better. 


SECTION  VIII 


CONCLUSIONS 


The  work  during  phase  II  of  the  project  has  accomplished  several  things  that  were  desired 
at  the  beginning,  though  some  results  were  not  anticipated: 

a .  An  Electronic  Shadowgraph  can  achieve  equal  or  better  accuracy  than  the  film  camera 
system  (using  geometric  methods),  provided  that  adequate  image  processing  be  accomplished 
and  the  centroid  and  second  moment  of  inertia  algorithms  be  used.  The  centroid  and  second 
moment  of  inertia  data  can  then  be  used  to  find  the  location  and  orientation  of  the  projectile  in  3- 
space  as  described  in  Section  VI. 

b .  Solid  state  cameras,  as  opposed  to  vidicon  cameras,  are  required  for  their  sensitivity 
to  work  with  the  existing  Eglin  spark-gap  as  part  of  an  electronic  shadowgraph.  Vidicon 
cameras  have  several  deficiencies,  the  worst  of  which  is  inadequate  sensitivity  for  the  Eglin 
application. 

c .  Solid  state  cameras  with  a  required  resolution  to  match  the  film  cameras  are  presently 
too  expensive  for  the  Eglin  Range,  and  are  forecast  to  remain  so  for  at  least  several  years. 

Cameras  with  a  2048  x  1024  resolution  could  cost  $50,000  each  in  100  quantities,  prohibitive 
for  Eglin . 

d .  Using  the  centroid  and  second  moment  of  inertia  algorithms  along  with  image 
processing  that  does  not  degrade  the  measurement  information  in  the  image,  solid  state  cameras 
of  1024  x  512  resolution  can  be  used  to  replace  or  augment  the  film  camera  system.  It  is 
strongly  recommended  that  the  electronic  shadowgraph  system  be  used  to  augment  the  film 
cameras,  and  that  they  both  be  used  in  tests  for  at  least  a  year  before  reducing  dependency  on  the 
film  cameras.  This  will  allow  collection  of  sufficient  data  to  compare  the  performance  of  the  two 
systems  and  to  attain  confidence  in  the  electronic  camera  system. 

e .  Solid  state  cameras  with  1024  x  5 12  resolution,  adequate  sensitivity  for  the  Eglin 
spark-gap,  including  an  A/D  converter  and  image  buffering,  and  communications  logic,  has  been 
estimated  to  cost  about  $7300  per  camera  (in  100  quantity)  by  a  leading  camera  manufacturer  and 
will  be  available  sometime  in  1988.  This  is  comfortably  within  the  stated  Eglin  goal  of  less  than 
$10,000  per  camera. 

f .  All  measurements  made  during  the  proof-of-principle  installation  at  Eglin  were  of  a 
coarse  nature  and  should  be  verified,  but  look  good  in  accuracy  comparisons  with  the  film 
system.  They  do  not  show  what  the  ultimate  accuracy  of  the  electronic  camera  system  could  be. 

g .  Significant  development  of  processes  to  adequately  prepare  our  images  for  the 
application  of  the  measurement  algorithms  must  be  done  before  higher  accuracy  measurements 
can  be  a  reality.  This  could  be  a  major  undertaking  but  may  result  in  some  new  developments  that 
could  be  applicable  to  several  Air  Force  and  DoD  test  and  measurement  problems. 
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h  .  No  detailed  study  of  how  to  locate  the  spin  pin(s)  usually  placed  on  projectiles  has 
been  done.  A  cursory  calculation  has  shown  that  a  2  mm  diameter  x  4  mm  long  pin  can 
adequately  be  located  with  a  1024  x  512  resolution  camera,  but  this  measurement  capability  must 
be  verified.  The  possibility  exists  that  locating  the  centroid  and  second  moment  of  inertia  of  the 
spin  pin  could  also  result  in  more  accurate  location  of  the  pin,  along  with  proper  three- 
dimensional  geometry. 

i .  Handling  of  winged  or  finned  objects  has  not  been  worked  out.  This  could  be  a 
problem  that  may  be  solved  by  proper  editing  of  the  image  using  a  priori  information  about  the 
model.  However,  much  work  and  experimenting  must  be  done  to  accomplish  it. 
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SECTION  IX 


RECOMMENDATIONS 


a .  It  is  feasible  to  begin  work  aimed  at  instrumenting  the  entire  Aeroballistics  Research 
Facility  with  electronic  shadowgraph  cameras  and  suitable  control  and  data  processing. 

b .  Continued  study  of  image  processing  techniques  for  measurement  applications  at 
Eglin  should  be  able  to  solve  those  related  problems.  High  accuracy  measurements  will  never  be 
achieved  until  suitable  image  processing  for  extracting  objects  from  images  can  be  done. 

c .  If  development  of  electronic  shadowgraph  range  instrumentation  is  to  continue, 
selection  of  a  solid  state  camera  and  associated  electronics  should  be  a  high  priority,  so  that 
camera  testing  and  characterization  can  be  done  prior  to  delivery  of  the  full  complement  of 
cameras,  or  perhaps  even  before  ordering  the  final  cameras. 

d .  Selection  and  purchase  of  an  image  processing  system/computer  should  be  completed 
very  early  to  allow  all  future  software  development  to  be  directly  applicable  and  useful  for  Eglin. 
Use  of  the  Vicom  system  and  the  vidicon  cameras  is  discouraged  because  of  its  numerous 
shortcomings  for  a  large  scale  installation,  not  the  least  of  which  is  its  clumsiness  with  camera 
interfacing.  That  system  is  also  over  four  years  old  and  has  demonstrated  problems  that  make  it 
unsuitable  for  good  software  development  and  long  term  use  and  support. 

e .  Completing  the  process  of  obtaining  measurements  with  the  prototype  system  and 
obtaining  final  location  and  orientation  results  using  Section  VI  equations  should  be 
accomplished  soon.  The  software  to  implement  the  equations  has  been  written  but  not  checked 
out.  This  is  an  important  step,  in  that  direct  comparisons  of  film  camera  and  electronic  camera 
results  can  be  made  for  many  shots  in  the  range.  This  will  allow  compiling  confirming  data  for 
the  electronic  cameras  and  methods,  and  will  be  required  for  proving  new  image  processing 
techniques  on  the  Eglin  images. 
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